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Abstract 


The effects of leading edge flaps on the aerodynamic characteristics of a low aspect- 
ratio delta wing are studied theoretically. As an extension of the classical cross-flow 
plane analysis and in order to include separated shear layers, an analogy between 
three dimensional steady conical and two dimensional unsteady self-similar flows 
is explored. This analogy provides a simple steady-unsteady relationship. The 
criteria for the validity of the steady-unsteady analogy are also examined. 

Two different theoretical techniques are used to represent the separated shear 
layers based on the steady-unsteady analogy, neglecting the trailing edge effect. 
In the first approach, each vortex system is represented by a pair of concentrated 
vortices connected to the separation points by straight feeding sheets. The analysis 
is carried out in the cross-flow plane by mapping the wing trace, by means of the 
Schwarz-Christoffel transformation, into the real axis of the transformed plane, 
allowing for exact satisfaction of boundary conditions. Although the model does 
not yield satisfactory comparison with experiments, the effects of leading edge flaps 
on delta wing aerodynamics are observed qualitatively. 

In the second approach, the vortex cloud method is adopted for simulating the 
flow field in the cross-flow plane. The separated shear layers are replaced with a 
cloud of discrete vortices and the boundary element method is employed to represent 
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the wing trace by a vorticity distribution. Combined with the steady-unsteady 
analogy, the vortex cloud method constitutes a numerically efficient methodology 
with good definition of separated shear layers and good prediction of surface pressure 
distribution. A simple merging scheme is used to model the core region of the 
vortical flow as a single vortex by imposing a restriction on the shear layer rotation 
angle. The results are compared with experiments and with results from 3-D panel 
calculations. 

In agreement with experimental evidence, it is found that flap deflection causes 
a decrease of lift and drag due to partial suppression of the leading edge vortical 
flow, while the lift-to-drag ratio increases due to the thrust component of the force 
acting on the flap surfaces. 


Acknowledgment 


This research was supported by NCC 2-55. 



Nomenclature 


a(x ) total half span length 

a,j geometric influence coefficient 

c speed of sound 

D drag 

e unit vector 

e total energy 

F complex potential in the cross-flow plane 
/ wing cross-section shape function 
G normalized vorticity distribution function inside core 
H total head 

i internal energy, \/— 1 

k ratio of main wing span length to total span length 
L lift 

lj length of j vorticity panel 

M Mach number 

m similarity coefficients 

N normal force 

p static pressure 

R resultant force, characteristic length scale 

r radial direction in cross-flow plane 
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r conical coordinate, Prandtl’s similarity coordinate 
S wing and vortex sheet geometry 

T characteristic time scale 

t time coordinate 

U chordwise component of free-stream velocity 
u velocity 

u |u| 

v disturbance velocity component 

v complex conjugate of the velocity in the cross-flow plane 
x vector position 

x chordwise coordinate of delta wing 

x,y,z stretched wing reference coordinate system through 
Prandtl-Glauert transformation 
y, z rectangular form of conical coordinates 
a angle of attack 

fa inside angle at the i th vertex of the polygon 
X vorticity distribution function inside core 
8 flap deflection angle, dirac delta function 

8{j kronecker delta function 

e half apex angle of main delta wing 

rj spanwise component of bound vorticity 

r vortex strength 


VI 
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specific heat ratio 

polygon i th vertex external angle 

viscous coefficient 


Ki 

V 

Q, vorticity on wing surface 

u vorticity 

u> | cZ? | 

$ total velocity potential 

<f> disturbance velocity potential 

<p dimensionless form of disturbance velocity potential 

\f f cross flow stream function 

p density 

q distance normalized with vortex core radius 
a a = y + iz, complex representation of cross-flow plane 
C vortex core radius 

0 0 = £ + i'd , complex representation of transformed plane 

02 , © 3 , © 5 , ©6 transformation coefficients 

f viscous stress tensor 

6 angular coordinate in cross-flow plane 

Subscripts 

c conical flow, core 
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F flap 



1 

h 

l 

l 

hinge vortex 

leading edge vortex 

P 

potential 

r 

rotational 

s 

separation point, body surface 

u 

2-D unsteady self-similar flow 

V 

vortex position 

w 

main wing 

x,y,z 

wing reference coordinate system 

OO 

freest ream 

( ) indicates complex conjugate 

O' 

, 

length scale normalized with a(x ) 

1 

’ 
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Chapter 1 
Introduction 


1.1 General Remarks on the Flow Field around Flapped Delta Wings 

The flow field on the leeside of low aspect ratio delta wings at incidence is char- 
acterized by a vortical flow originating at the leading edges, as shown in Fig. 1. 
This vortical flow is very stable up to a certain angle of attack and causes an en- 
hancement of lift, called vortex lift, by inducing low pressure regions on the upper 
surface. Generating a stable vortical flow above the wing surface is an efficient 
way of obtaining additional lift, 1 especially during maneuvering and in take-off and 
landing. At high angle of attack, however, the vortical flow also causes the following 
problems: 

• Large lift-dependent drag due to loss of leading-edge suction 

• Side force due to asymmetric formation of vortical flows 

• Pitch-up due to vortex breakdown at high angle of attack 

Fig. 2 shows typical concepts for handling these problems, as well as the basic 
functions of the concepts. 
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Leading edge flaps have been applied to the design of subsonic-transonic tactical 
delta wing aircraft, and to aircraft that require both supersonic cruise efficiency 
and high maneuverability. 2 On the other hand, the flow field about a wing with 
leading edge flaps is more complex than that of a wing without flaps due to the 
appearance of an additional vortical system originating at the flap hinge, as shown 
in Fig. 3. It has been shown that the effect of these hinge vortices is negligible, 
except at small angle of attack with large flap deflection angle. ® Various shapes 
of leading edge flaps have been studied experimentally in the past, in an attempt 
to optimize their influence on the aerodynamics of delta wings. 4-14 One of the 
effects of leading edge flaps is to increase the lift-to-drag ratio by reducing the 
wing drag. This augmentation occurs due to a thrust component of force acting 
on the flap surfaces. Fig. 4 shows the variation of lift-to-drag ratio in cruise for 
three different aircraft configurations: the F-16 aircraft, a precramp wing* aircraft 
of supersonic design, and a flapped delta wing, for a wide range of Mach number. 
Leading edge flaps also have the detrimental effect of partially suppressing the 
vortical flow, thereby reducing the total lift. These two phenomena have been 
observed experimentally'*’ 12-1 * 5 and theoretically. 17,18 

1.2 Previous Studies 

In view of the important influence of vortical flow on wing design, numerous studies 
have been conducted on the leading edge vortices around low aspect ratio delta 
wings. 20-51 These studies were aimed at understanding the fundamental physics 
of rolled-up shear layers and the effect of the vortical flows on the aerodynamic 
characteristics of the wing. Main topics of such studies were: 


* a cranked wing planform with twist and camber 



1. Introduction 


3 


• Mathematical modeling of rolled-up shear layers 

• Flow field analysis of the core regions of vortical flows 

• Three-dimensional (3-D) boundary layer calculations and prediction of sepa- 
ration lines 

• Vortex breakdown 

The analysis of the core region of the vortical flow is essentially a part of the study 
of the rolled-up shear layer. The main difficulties in the study of the core region 
arise from the mathematical treatment of the tightly coiled shear layers and the 
effect of viscosity. ^* 4- ^ The calculation of the 3-D boundary layer^®^® and the 
prediction of separation lines 51 ’^ 2 can not be made until the flow field external 
to the boundary layer, a flow field composed of rolled-up shear layers, has been 
determined. A description of the core regions combined with trailing edge effects is 
also necessary for an improved understanding of vortex breakdown. In this regard 
very little progress has been made in the analysis of this problem, however. 

The first successful model for the vortical flow was a straight-feeding-sheet model 
used by Brown and Michael, as shown in Fig. 5. 2 ®’ 2 * Other attempts have used ana- 
lytical approaches, 2 ® -2 ® vortex cloud methods,®® - ®® and 3-D panel methods. 4 ® -42 
More recently, with the aid of high computing power, simulations of separated vor- 
tical flows have been carried out using numerical solutions of either the Euler or 
Navier-Stokes equations. 4 ^ - ®® 

Studies on control of vortical flow were started in late 1970’s, after some un- 
derstanding of the physics of vortical flow had been achieved. At the beginning, 
attempts to suppress the vortical flow were made when it was recognized that con- 
ventional control surfaces could not operate well in this type of flow. Leading 
edge flaps were found to be quite effective in providing vortical flow control, as was 
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revealed in experiments conducted with complex-geometry delta-wing aircraft. 1 ^’ 1 ® 
Those experimental results showed that leading edge flap deflection was effective in 
promoting an attached flow configuration by suppressing the formation of leading 
edge vortices. The tabbed-vortex-flap* has also been considered to avoid asym- 
metric vortex breakdown with sideslip angle. Both plain leading edge flaps and 
tabbed-vortex-flaps caused significant increments in lift-to-drag ratio. 

In the experimental study by Rao et various shapes of leading edge devices 

were suggested to reduce lift-dependent drag. Frink 19 performed a study of the 
effects of leading edge flaps by using a 3-D free-vortex-sheet method.** Hoffler et 
al . 14 investigated the aerodynamic influence of tabbed-vortex-flaps and plain lead- 
ing edge flaps, both experimentally and using a 3-D panel method. Marchman 11,12 
performed similar experiments for different delta wing configurations in an attempt 
to confine the vortical flow to the vicinity of flap surfaces in order to reduce drag, 
while keeping the lift almost constant. The effect of a vortex fence, which is an 
alternative form of leading edge flap, as shown in Fig. 2, has been investigated by 
Reddy using a 3-D Free Vortex Sheet Method. 19 Lee 3 made experiments with a 
flapped delta wing with conical geometry, including flow visualization, to observe 
the influence of the hinge vortices. 

1.3 Present Study 

In this study, two different theoretical techniques are used to assess the effects of 
leading edge flaps on the aerodynamic characteristics of slender delta wings. The 

* An alternative form of leading edge flaps; a second leading edge flap is attached to the main 
flap with different deflection angle. 

** Frink’s 3-D panel calculations made use of only one vortical system, existence of the hinge 
and trailing edge vortex systems being neglected. 
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first technique uses the simplest vortex model in a plane normal to the wing surface 
(cross-flow plane) to explore the general trends of the aerodynamic coefficients of 
flapped delta wings. The second technique improves on the first by incorporating 
an in viscid vortex cloud simulation of the shear layers. Both approaches ignore the 
presence of the trailing edge. 

The foundation of both approaches is based on the mathematical analogy be- 
tween conical steady and two-dimensional (2-D) unsteady self-similar flows, which 
results in a simple relationship between time and space coordinates. With this anal- 
ogy (steady-unsteady analogy) a steady conical flow field can be reduced to a 2-D 
unsteady self-similar flow in the cross-flow plane. This steady-unsteady analogy 
under the conical assumptions yields more accurate results in the case of supersonic 
flow,^ 1 ^ since the whole flow field is conical, whereas for subsonic flow the analogy 
ceases to be valid near the trailing edge and the apex. The maximum flap deflection 
angle for which the steady-unsteady methodology remains valid is found by using 
a simple geometrical relationship. A deflection beyond that limit would cause the 
leading edge vortex system to move under the wing surface, a situation that cannot 
be accounted for in the context of a cross-flow plane analysis. For both techniques, 
the effect on the lift-to-drag ratio of the resultant forces acting on wing and flap 
surfaces is derived to first order using a coordinate transformation. 

The two approaches are now described in more detail. The first approach consti- 
tutes an extension of the vortex model first suggested by Brown and Michael. The 
study used here simulates each vortex system by a pair of concentrated vortices 
with straight feeding sheets. The unknowns in the problem, vortex strength and 
position, are found from force equilibrium conditions imposed on the singularity 
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systems. The Schwarz-Christoffel transformation is used to allow for exact sat- 
isfaction of boundary conditions. Although the model does not give satisfactory 
results on vortex position and pressure distribution the global quantities, such as 
lift, drag and normal force, are satisfactorily calculated. Particular attention is 
paid to the following aspects: 

• Influence on lift-to-drag ratio of forward components of the forces acting on 
the flap surfaces 

• Variation of vortex position and strength with flap deflection angle 

• Contribution of vortical and attached flows* to the aerodynamic forces 

In the second approach, the vortex cloud method, the steady-unsteady analogy 
plays a more explicit role. Assuming that geometric and aerodynamic slenderness 
are satisfied, this analogy reduces a conical steady flow field to a 2-D, unsteady self- 
similar flow through the relationship x = Ut, where x is the chordwise coordinate 
and U is the chordwise component of the free-stream velocity. Using this analogy 
the vortex cloud method, an intrinsically 2-D approach originally developed to de- 
scribe separated flow past bluff bodies, 55-76 is modified to simulate a conical flow. 
This approach leads to a numerically efficient methodology with good definition 
of shear layers. Another motivation for the use of this approach is the parallel 
between the physics of planar shear layers and the delta wing rolled-up shear layer, 
as suggested in Fig. 6. 77 

Rolled-up shear layers are replaced by discrete vortices having finite core size, 
so that all of the rotational and viscous effects are assumed to take place within 
the core, while outside the core the flow field is governed by conventional potential 

* As will be discussed later, the lift of a delta wing with separated flow may be expressed as the 
sum of a contribution that would have been caused by a wing with attached flow and a component 
due to the separated flow. 
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theory. Viscous diffusion of each vortex is taken into account through a formulation 
derived by Prandtl. 78 In conical flow the “vortex jungle ”, which usually forms in 
the wake region in 2-D flow simulations, would represent the region of tightly coiled 
shear layers, or core region of the delta wing flow field. Such a representation of the 
core region as a “vortex jungle” is inappropriate in conical flow due to the appear- 
ance of the high axial velocity, which would violate the aerodynamic slenderness 
assumption and could not be modeled in a 2-D flow context. 7 ^ To improve the 
computational efficiency as well as to avoid this unphysical core representation, a 
merging scheme between neighboring vortices and near the core region is employed. 
A quasi-steady Kutta condition is imposed at leading edges and flap hinges by fix- 
ing the separation angle of new vortices generated at each computational step. If 
the Kutta condition were unsteady the new vortex strength and position would be 
coupled with the strengths and positions of other vortices in the shear layers and 
those of vortices representing the wing surface, leading to a matrix inversion oper- 
ation for computing the bound vortex strengths. With this treatment of the 

Kutta condition matrix inversion, which requires a great computational effort, can 
be avoided at each iteration by decoupling the Kutta condition from the system of 
equations. 

In the context of the vortex cloud method, attention is focused on the following 
aspects: 

• Distortions of rolled-up shear layer with flap deflection angle 

• Pressure distribution 

• Variation of vortex position and strength with flap deflection 

• General features of trailing edge wake structure 

The last item is explored under the assumption that the conicality of the flow field 
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is sustained up to the trailing edge. These results are compared with experimental 
and other theoretical studies, including results from 3-D panel computations. 


1.4 Overview 

This work is divided into six chapters. Chapter 2 contains the derivation of the 
basic mathematical equations, including the analogy between conical steady and 
2-D unsteady self-similar flow, the range of flap deflection angle where the steady- 
unsteady analogy is valid, and the lift-to-drag ratio of a flapped delta wing to the 
first order. Chapter 3 describes the first approach together with the fundamental 
concepts of the Brown-Michael vortex model and the Schwarz-Christoffel transfor- 
mation. In Chapter 4 the basic mathematical background of the vortex cloud 
method for a 2-D unsteady flow is derived and adapted to conical flows, and nu- 
merical implementation is carried out. At the end of these chapters the results 
of each method are presented and the mathematical and numerical limitations of 
each technique are discussed. Chapter 5 shows the comparison of the results from 
the two methods with those from other analyses. Chapter 6 contains discussions 
on the relative accuracy of each of the techniques implemented here as compared 
with experiments and other methods of calculations, as well as recommendations 
for future work. 




Chapter 2 
Basic Mathematical Analysis 


The basic mathematical analysis to treat steady 3-D conical flows is derived. The 
linearized potential equation is derived from the Navier-Stokes equations, and trans- 
formed through the Prandtl-Glauert transformation. The assumption of conicality 
of the flow field around delta wings is examined for the subsonic and supersonic 
cases. Criteria for applicability of the cross-flow plane analysis are obtained through 
the mathematical analogy between a steady 3-D conical flow, which is governed by 
the linearized potential equation, and a 2-D unsteady self-similar ideal flow. The 
range of validity of the steady-unsteady analysis is examined. The fundamental 
equations for the two approaches, the straight-feeding-sheet and vortex cloud mod- 
els, are presented on the basis of the analogy. An expression for the lift-to-drag 
ratio of flapped delta wings is derived to first order. 

2.1 Linearized Potential Equation 

As described in Chapter 1, the flow field around a delta wing is composed of the 
boundary layer on the wing surface, the rolled-up shear layers originating at the 
leading edge, the secondary vortex systems originating in flow separation from the 
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upper surface, and the trailing edge wake, as shown in Fig. 1. The equations 
governing this flow field are: 


Conservation of mass: 

Conservation of momentum: 

d(pu) 


dp 

dt 


+ V • (pu) = 0 


dt 


+ V • (puu) = — Vp + V • f 


Conservation of energy: 
dpt 


+ V-(/>eu) = -V-(pu) + V-(f-u) 


( 2 . 1 . 1 ) 


( 2 . 1 . 2 ) 


(2.1.3) 


where f is the stress tensor which, using Stokes hypothesis for the molecular vis- 
cosity, can be expressed as 


( dui duk 2 du\ 

Tik ~ d ^ + ~dTi ~ 3 6ik dTi 

and e is the total energy per unit mass 

e = i + u 2 /2 

The boundary conditions are 

u = 0 at S(x, y,z) = 0 
u = U as x — * oo 


(2.1.4) 

(2.1.5) 

(2.1.6a) 

( 2 . 1 . 66 ) 


where S ( x,y,z ) = 0 represents the wing surface geometry. 

Usually an aircraft operates in the range of high Reynolds number, where the 
thickness of the boundary layer and the shear layer become so small that all the vis- 
cous effects are confined to thin regions. Thus, the viscous effects can be neglected 
in most of the flow field. The boundary condition on the wing surface is 


u • n = 0 


(2.1.7) 
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where n is the outward normal vector on the wing surface. 

In the absence of viscous effects the rolled-up shear layers are infinitely thin with 
boundary conditions requiring that there must be no pressure jump across them. 
These concepts lead to an inviscid model for the flow field governed by the Euler 
equations. The flow field is still rotational, in which case the velocity can be 
decomposed into two components as follows 

u = Up + u r (2.1.8) 

where u p and u r are the potential and rotational part of the velocity respectively. 
These two components have the following properties 

u p = V$ or V x Up = 0 (2.1.9a) 

V x u r = w (2.1.96) 

where $ is a scalar potential and Cj is the vorticity. 

Neglecting the viscous terms Eqs. (2.1.1) and (2.1.2) become 

7^7 = ~/>V 2 $ - pV • u r (2.1.10) 


V 


(?*, [P dp) 

[dt + 2 + J Pao p) 


dn r 

~dt 


U p 

• U r + — 


“l - Up) X UJ 


( 2 . 1 . 11 ) 


Assuming that the flow field is homentropic the following relationships hold 

D [P dp _ d / fP dp\ Dp 1 dp Dp 

Dt Jpoo P d P \Jpo o p ) Dt p dp Dt 

_ c 2 Dp 
p Dt 


-C 2 (v 2 $ + V • Ur) 


( 2 . 1 . 12 ) 


2. Basic Mathematical Analysis 
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Combining Eqs. (2.1.11) and (2.1.12) there results 
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v2 *-?@ + s' V4|2 + V4 - v (5 iV412 )) 

= —V • u r 

•V^-£ + A -dx + (V$ + u r )-A^ 

(2.1.13) 

where 

(2 \ 

V$-u r + yj + (V$ + u r )xw (2.1.14) 

The left-hand side of Eq. (2.1.13) is exactly the same as the full potential equation 
for a compressible gas. The right-hand side can be interpreted as the source of 
rotational effects. 

Let 

$ = JJx + U tana y + <j) (2.1.15) 

where a is the angle between the free-stream and the x direction and 

U r = (urn ^ryi ^rz) (2.1.16) 



Assuming the flow field is steady and all the disturbance velocities are small com- 
pared to the uniform free-stream velocity Eq. (2.1.13) gives, after neglecting terms 
of second and higher order in a, ^ and u r 

d 2 <f> 1 / d 2 (f> d 2 <f>\ 

dx 2 1 - Moo 2 \d? dz 2 ) 

du rx 1 f du r y du rz \ 

dx 1 - Moo 2 V dy dz ) 


(2.1.17) 
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where M 0 o is the free-stream Mach number. 


Applying the Prandtl-Glauert transformation for subsonic flow 

x = x, y = y yjl- Moo 2 , z = z yj 1 - Moo 2 (2.1.18) 

to Eq. (2.1.17) there results 

d 2 <^> 5V , du rx 1 f du ry durz\ ( . 

dx 2 3y 2 dz 2 dx \/l - Moo 2 V ^ dz / 

which is a form of the Poisson equation. 


As mentioned before, all of the rotational effects are confined to thin layers. Green’s 
theorem will be applied by treating the shear layers as part of the boundaries. 
Green’s theorem states that a potential function can be obtained through a suitable 
distribution of singularities on the boundaries, 80-82 in this case the wing surface 
and the shear layers, as follows 


^ = ib/R v ^" d5 - ^M n ' v R ) d5 < 21 - 2 °) 

where S represents the surface of the boundaries and n is the outward normal 
vector. 


Eq. (2.1.20) shows that the velocity potential <f> is due to a surface distribution of 
sources and doublets. The density of sources is — V<f>- n per unit area. The density 
of the doublets is <f> per unit area and the doublet axis is directed along the inward 
normal to the surface. 


Once the boundaries are replaced with a suitable distribution of singularities, the 
governing equation of the flow field outside the boundaries is 

d 2 (f> d 2 <j> d 2 6 

aF + W + aW ~ 0 (2 ' L21) 

For supersonic flow the Prandtl-Glauert transformation is 


x = x, 


y = y \/ Moo 2 - 1, 5 = z \J Moo 2 - 1 


( 2 . 1 . 22 ) 
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and the governing equation becomes 

d^4> d^4> du rx 

dx 2 + dxp 1 dz 2 dx 


1 ( ® u ry 

VMoc 2 - 1 v dy 



(2.1.23) 


In the supersonic case the confined rotational effects can also be replaced with 


suitable singularities,®*’®'* leading to a statement similar to Eq. (2.1.20). The 


governing equation outside the singular regions is in this case 

_ - n 
dx? dy 2 5z 2 


(2.1.24) 


Eqs. (2.1.21) and (2.1.24) show that the linearized potential equation can be 
reduced through the Prandtl-Glauert transformation to 


, a 2 ^ dH a 2 ,» n 

± 5F + W aF “ 0 


(2.1.25) 


which implies that the governing equation in case of subsonic flow is transformed to 
that of an incompressible flow and in supersonic flow to that of another flow with 
free-stream Mach number equal to \/2. 


2.2 Cross-Flow Plane Analysis 
2.2.1 Delta Wing Flow Conicality 

The steady-unsteady analogy, which will be discussed later, is conducted under 
the assumption of conicality of the flow field. The conical assumption for a flow 
field requires that the wing geometry be conical and that there should not exist 
any characteristic length scales.®'* In a conical flow field, velocity and pressure are 
constant along rays emanating from the wing vertex. 

Conicality of the flow field around a delta wing is better discussed using the small 
disturbance potential equation expressed in physical space. If the shear layers are 
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treated as a part of the boundaries, Eq. (2.1.17) becomes 

I? + T^J iw + £*) = 0 (2 ' 21) 

For a supersonic flow the value 1 — Moo 2 becomes negative, resulting in a wave 
equation, with solutions propagating in the positive x direction. Hence, there is 
no effect of the length scale given by the trailing edge, since all of the disturbances 
as a result of the wing propagate only in the downstream direction. 


For a finite delta wing subsonic conicality is an approximation that stems from 
ignoring the presence of the trailing edge and the singular nature of the apex, and 
must be interpreted as a local solution near the wing apex region. Hence, conicality 
for subsonic flows means that the area of interest is removed from the trailing 
edge and at some distance from the apex. This is consistent with experimental 
observation.'* 


2.2.2 Analogy between 3-D Conical and 2-D Unsteady Self-Similar Flows 

In this section, as a extension of the classical slender body theory, 53,54,81,82 which 
does not consider separated shear layers, a comparison is made between a three- 
dimensional, steady flow field governed by Eq. (2.1.25) and growing conically in 
space, and a two-dimensional, self-similar ideal flow field growing in time, as illus- 
trated in Fig. 7. 

For a conical flow field around a delta wing with incidence angle a , the equation 
expressing the velocity potential can be written 

$ c (x,r,0) = x U (1 + rtana sin# + <£> c (r,0)) (2.2.2) 

where r = r/x, and <p c is the dimensionless form of the disturbance velocity 
potential defined as 

n o \ 
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The geometry of the vortex sheet and the wing surface can be also written 


S(x,r,9 ) = r - x / c (0) = 0 


(2.2.4) 


where x f c represents the cross-section of the wing and the rolled-up shear layer 
growing linearly in the x direction. 

In terms of <£> c , r and 9 , Eq. (2.1.25) becomes 

_2 5 2 v? c 1 d 2 <p c 


1<Ve + (i ±1*)^ + a 


= 0 


(2.2.5) 


r dr ' v ’ dr 2 ' 89 ^ 

where + and — correspond to the subsonic and supersonic case respectively. 

The condition for the boundaries, the wing surface and the vortex sheet to be stream 
surfaces is 

V$ c • VS = 0 (2.2.6) 

Combining Eq. (2.2.6) with Eqs. (2.2.2) and (2.2.4) gives 




. n dfc 
tan a sin 9 + — 
or 


tl 

fc 


I J c I ± q | 1 

+ — tana cos0 + 

r o9 


(2.2.7) 


= 0 


where V 3 ./U = <p c ~ r(dip c /9r) 

The boundary condition that there be no pressure difference across the vortex sheet 
can be formulated by invoking Bernoulli’s equation 


P-Poo 

(pvy 2 ) 


^1 + tan 2 a^ — ^1 + — (tana sin0+-^ 


— (tan a cos 9 + 


U £\ 2 

u) 


( 2 . 2 . 8 ) 


where p is the pressure in the transformed space. 
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In terms of <p c , r and # Bernoulli’s equation can be written as 



(2.2.9) 


The value of this pressure must be the same on either side of the vortex sheet. 

For an unsteady, self-similar, 2-D flow, Prandtl’s analysis 22 showed that the 
vortex sheet has a shape that belongs to a self-similar family. Following his analysis 
with the assumption that the free-stream velocity causing separation is U tana, the 
velocity potential and the geometry of the boundaries, including the vortex sheets, 
can be written 

d 2 / ^ \ 2m- 1 

$ u (t,r,0) = Ur tana sin0 + — f — J <fu(r,0) (2.2.10) 


S(t,r,0) = r - R (!) /»(«) = 0 


( 2 . 2 . 11 ) 


where 



( 2 . 2 . 12 ) 


and is the dimensionless form of the disturbance velocity potential defined as 

<p u = H T (2.2.13) 

v u R 2 (t \2m-l v ’ 

t yr) 


and T and R are constant reference values for time and length, m is a parameter 
which designates a particular member of the family of the vortex sheets. f u is 
the collapsed geometry of the family of the vortex sheets and the wing geometry 
when those are represented in terms of the similarity variable r. In this study the 
geometry of the wing surface is also growing with time. 
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For the 2-D unsteady case Laplace’s equation can be written in terms of ip u as 

ldv? u d 2 (pu , 1 d 2 Vu 


+ 


r dr ' dr 2 r 2 dQ 2 
The kinematic boundary conditions of the vortex sheet and wing surface are 8 ^ 

dS 


= 0 


(2.2.14) 


dt 


+ VS.-VS = 0 


§ + V **- V5 = 0 


where subscripts i and e represent each side of the vortex sheet, and 


V$ t - = 


V*| s + 2 




7 

2 


(2.2.15a) 

(2.2.156) 


(2.2.16a) 

(2.2.166) 


where V$| s is the value on the vortex sheet, and 7 is the vorticity intensity of the 
sheet. 

Combining Eqs. (2.2.15) and (2.2.16) there results 

~ + V$| 5 -V5 = 0 (2.2.17) 

In terms of r and <~p u this equation becomes 


m f u - 


+ 



fx i 

fu 


0 


(2.2.18) 


The dynamic boundary condition on the vortex sheet is that the pressure must 
be the same on either side of the sheet. The pressure can be obtained from the 
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unsteady Bernoulli’s equation 
P ~ Poo 


-C p - 


(pU 2 ) / 2 


t 2 , 1 2 2 
= -tan a + 2 -^— + u r - vq 

ot 


(2.2.19a) 



All the governing equations for the two cases, Eqs. (2.2.5), (2.2.7), (2.2.9) (the 
steady case) and Eqs. (2.2.14), (2.2.18), (2.2.19b), (the unsteady 2-D case), are 
written in nondimensional form. The unknowns become the scalar functions y? c 
and <p u , the collapsed form of the boundary geometry, including that of the vortex 
sheet, f c and f u for both cases. If the governing equations and auxiliary conditions 
were identical, the solutions would be the same in both cases. The conditions under 
which such is the case provide the criteria for a mathematical analogy between a 
steady conical flow and an unsteady, self-similar, 2-D flow. The mathematical 
requirements for such analogy to exist are as follows: 

Laplace’s equation for the two cases, Eqs. (2.2.5) and (2.2.14), becomes formally 
identical if r = 0 . Hence the requirement 


r < 1 


( 2 . 2 . 20 ) 


The form of the boundary conditions for the vortex sheet and the wing surface, Eqs. 
(2.2.7) and (2.2.18), become the same if v x /U = 0,m = 1 and R/T = U . Hence 
the requirements 
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m = 1 (2.2.22) 

f = v (2.2.23) 

Conditions (2.2.22) and (2.2.23) also guarantee that Bernoulli’s equation on the 
vortex sheet, Eqs. (2.2.9) and (2.2.19) becomes identical. The relationship between 
x and t is given by 

x = Ut (2.2.24) 

The analogy shows that the governing equations for a 3-D conical flow reduce to 
those for a 2-D unsteady self-similar flow if following criteria are met 

• r <C 1 

• Vx <C U 

These conditions mean that the conical flow field must be geometrically and aero- 
dynamically slender. 


2.3 Equations for the Straight-Feeding-Sheet and Vortex Cloud Models 

2.3.1 Space-Time Transformation with Conical Flow Assumptions 

With the analogy described in the previous section, a 3-D conical flow field can 
be obtained from the results of a 2-D unsteady self-similar flow field. The two 
approaches, which will be described in the next two chapters, are formulated on 
the basis of solving a 2-D unsteady flow in the cross-flow plane instead of solving 
3-D conical flow directly. The solutions of the 2-D unsteady flow should reach a 
converged state when they are represented in terms of the similarity variable, due to 
the requirement of the self-similarity. From the relationship given by Eq. (2.2.24), 
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the time derivative can be rewritten as 


T = U 7~ 
at ax 


(2.3.1) 


The length scales in the cross-flow plane should be normalized with a characteristic 
length, which is proportional to t , or from the relationship x = Ut, proportional 
to x . In both approaches all the length scales are normalized with the total span 
length a(x). From Fig. 8, a(x ) is given by 

x tan e 


a(x ) = 


(2.3.2) 


where k is the ratio of the main wing span length to the total wing span length. 
Combined with the conical assumption, Eq. (2.3.1) reduces to 


U ^ = jj - U tanc S 

dx x ka(x) 

(2.3.3) 

U tane w 
k 

where E stands for any quantity, growing linearly with x, associated with the 
kinematics of the flow field. As a consequence of this, Y,' = Y/a(x) = constant. 


2.3.2 Basic Equations for the Straight-Feeding-Sheet Model 

In this approach the vortex structure is specified as consisting of a straight vortex 
sheet, emanating from the leading edge and flap hinges, and concentrated vortex 
cores, as shown in Fig. 9. The steady-unsteady analogy of the previous section 
showed that the analysis can be carried out in a time dependent fashion in the 
cross-flow plane, by solving Laplace’s equation within prescribed boundaries, Eqs. 
(2.2.14) and (2.2.17). The vortex sheets are part of the boundaries. In the 
straight-feeding-sheet model Eq. (2.2.17) cannot be satisfied. This means that 
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the singularities representing the separated vortex sheets cannot be made force- 
free at every point. Instead of Eq. (2.2.17), a global force-free requirement for 
the vortex+sheet system is introduced. By representing the cross-flow plane with 
complex notation, the force per unit length acting on the j th vortex+sheet system, 
as represented in the unsteady cross-flow plane formulation is 8 ^’ 88 

P iTj + P ~1I^~ ^ vj ~ = ( 2>3>4 ) 

where the cross-flow plane has been replaced by the complex representation a — 
y + iz. Tj is the vortex strength and a v j = y v j + iz v j is the vortex position, v v j 
is the velocity at the vortex position excluding its own induced velocity, and a s j is 
the origin of the vortex sheet. 

The first term is the force acting on the vortex core, the second is that acting 
on the vortex sheet. The term ( dcr v j/dt ) — v v j is the relative velocity causing a 
Joukowsky force acting on the concentrated vortex core. 87 Equilibrium requires 
that this force be balanced by forces acting on the sheet. This will be imposed on 
each vortex+sheet separately. Eq. (2.2.14) is rewritten in complex notation as: 

4 

V 2 F = iT j % - °j) ( 2 - 3 - 5 ) 

3 - 1 

where F is the sum of the complex potential due to the vertical component of 
the free stream velocity, concentrated vortex and the wing trace, and 6 (a — crj) 
represents the dirac delta function, required to represent a singular vortex. 

The basic equations for the straight-feeding-sheet model are the global force-free 
conditions: 


fj = 0 


j = 1,2,3, 4 


(2.3.6) 
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2.3.3 Basic Equations for the Vortex Cloud Model 

As justified by the steady-unsteady analogy of Section 2.2.2, the approach will 
be implemented in a 2-D unsteady form. Moreover, the vortex cloud method is 
a well-established technique to simulate a 2-D unsteady separated flows. In this 
approach the shear layers are replaced with a cloud of discrete vortices. The 
governing equation and the boundary condition given by Eqs. (2.2.14) and (2.2.17) 
were derived by treating the separated shear layers as a part of the boundary. In 
contrast, in the vortex cloud method the shear layers are not considered as a part of 
the boundary, instead, the shear layers are computed by tracing a large number of 
vortices, whose intensities and positions become known as the calculation proceeds. 
If the shear layers, represented by a vorticity distribution u>, are included in the 
field rather than in the boundaries, Eq. (2.2.14) takes the form 

V 2 # = -u (2.3.7) 


where ’3' is the stream function. 

The reason for the use of a stream function instead of a velocity potential is that, 
in terms of the stream function, the boundary value problem in the cross-flow plane 
is of Dirichlet type, and hence well-posed. 

The tangency condition represented by Eq. (2.2.17) is now replaced by the 
requirement that the vorticity representing the shear layers be convected with the 
local velocity 

Duo 

Di - ° < 2 ' 3 - 8 > 

which implies that the vorticity of each particle is not changed during the convection 
and represents the conservation of vorticity. 
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In addition, the total circulation in the flow field, including the bound vortex 
due to the wing trace and the vortices distributed on the shear layers, should be 
conserved: 


d(Tw + ?Sl) 
dt 


(2.3.9) 


where T and T gi represent the circulation bound to the wing trace and that of 
the shear layers, respectively: 

T w 

r SL 

where v s is the velocity on the wing surface, dl follows the contour around C along 
the wing trace, and S is the area including all the vortices representing the shear 
layers, as shown in Fig. 10. 


= l v ‘ 
-II 


uidS 


(2.3.10a) 

(2.3.106) 


The boundary conditions on the wing trace and the far field condition in terms of 
4' are 


^ = constant on wing surface (2.3.11a) 

\I r = $oo as r — ► oo (2.3.116) 

Details of implementation will be presented in Chapter 4. 

The fundamental equations for the vortex cloud model are then Eq. (2.3.7) with 
boundary conditions (2.3.11), and Eqs. (2.3.8) and (2.3.9). 


2.4 Validity Range of the Cross-Flow Plane Analysis 


2.4.1 Slenderness of Flow Field around Delta Wings 

In the flow field around delta wings having half apex angle e and at angle of attack 
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a , the disturbances in the cross-flow plane, such as the rolled-up shear layers, remain 
within a region whose radius is one to two times the half span length, except when 
the angle of attack is extremely high. Thus the region of interest in this study is 
limited by a radius given by 

r ~ x tan e<i (2.4.1) 

since for a low aspect ratio delta wing the half apex angle is small. This satisfies 
the condition mentioned in Section 2.2.2 

r < 1 (2.4.2) 

This also implies 

0(f) = 0(e) (2.4.3) 

Experimental observations 7 ^ have shown that the condition v % <C U is also satisfied 
except in the core region of the rolled-up shear layers, where the axial velocity is 
larger than free-stream velocity. Hence, it can be assumed that the flow field 
around a slender delta wing at moderate angle of attack satisfies the geometric and 
aerodynamic slenderness conditions. 

2.4.2 Relationship between a and e 

Although the steady-unsteady analogy developed in Section 2.2.2 does not impose 
any limitations on the relative magnitude of a and e, as long as both are small, 
experimental observations indicate that the ratio a/e must be within certain upper 
and lower bounds for the steady-unsteady analogy to hold. 

For large a/e the rolled-up shear layers do not form stable conical shapes but 
develop in a Karman- vortex-street pattern, thereby destroying conicality. If a 
is of order smaller than that of e, the shear layers roll up tightly within a small 
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region, thus becoming immersed in a region of dominant viscous effects. The 
viscous effect provides the characteristic length scale - \ffi/{pU ) and as a result the 
flow field cannot be conical, which requires no characteristic length scale. Thus, 
the steady-unsteady analogy is valid if a and e are of the same order. 

2.4.3 Limit of Flap Deflection Angle 

The cross-flow plane analysis in the framework of the steady-unsteady analogy, as 
described in Section 2.2.2, implies that the significant velocity component causing 
separation along leading edges is the vertical component of the free-stream velocity 
in the cross-flow plane. The maximum flap deflection for which such analysis is 
valid is related to the velocity component normal to the flap surface. This can 
be shown by considering that, if that component were directed towards the upper 
flap surface, as shown in Fig. 11, the leading edge vorticity system would tend to 
be located under the wing. The cross-flow plane analysis cannot, by definition, 
exactly account for the actual orientation of the velocity component normal to the 
flap surface. Since the discrepancy between the actual velocity component normal 
to the flap surface and the one assumed by the cross-flow analysis increases with 
flap deflection, this method is expected to fail for flap deflection beyond a certain 
range. An assessment of that range can be made by considering the condition 
under which the velocity component normal to the flap surface vanishes (a situation 
unaccountable for in the cross-flow analysis.) This assessment is facilitated if the 
velocity component normal to the flap surface is assumed to be due to the free- 
stream velocity alone, that is, if the perturbation velocity of the wing is neglected. 
This is plausible if the angle of attack and apex angles are both small. 

The outward normal vector, n, on the flap surfaces can be obtained through 
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the coordinate transformations as shown in Fig. 12. e x ,ey,e 2 denote 
vectors along axis x,y,z. The relationships between the unit vectors (e 
and (e Xl ,e yi ,e 2l ), and between (e Xl ,e yi ,e 2l ) and (e X2 , e y2 , e 22 ) are 

( e Xl \ / cos e 0 sin e \ 

e !h = ? i 0 

e 2l / \-sm e 0 cos e/ 



( e * 2 \ / 1 0 0 \ 

I e y 2 J = 10 cos 8 sin 8 I 

\e Z2 J \ 0 — sin 8 cos 8 ) 

Identifying n with the outward normal vector on the upper flap surface, 
seen from Fig. 12 that 



n = e y2 


Then 

n = cos 8 e yi + sin 8 e 2l 

= — sin e sin 8 e x + cos 8 e y + sin 8 cos e e 2 


The free-stream velocity has two components such that 


V = U e x + U tan a e y 

The component of the free-stream normal to the flap surface is then 
U n = U • n = -U sin 8 sin e + U tan a cos 8 


The condition that U n > 0 is 


tan a c 

— > tan 8 

sin e 


the unit 

E 5 Gy, e 2 ) 

(2.4.4) 

(2.4.5) 
it can be 

(2.4.6) 

(2.4.7) 

(2.4.8) 

(2.4.9) 
(2.4.10) 


The boundary shown in Fig. 13 represents the limit where the component of the 
free-stream velocity normal to the flap surface vanishes for e = 15 deg. 
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In this section, a general formulation for the lift-to-drag ratio of a flapped delta wing 
is derived. Experiments 3,4,15,16 show significant increment of lift-to-drag ratio with 
flap deflection. The rotation of the leading edge flap about the flap hinge axis causes 
the normal vector to have a component opposite to the flight direction, as derived 
in the previous section. The first term in the last line in Eq. (2.4.7) indicates 
the forward component of the normal vector. The associated force component 
constitutes the thrust force and is of higher order than the other two, under the 
assumption that e is small. The fact that flap deflection causes suppression of the 
vortical flow 6 can be interpreted as the vortical flow being confined to the vicinity 
of the leading edge flaps. Hence, most of the pressure peak region induced by the 
vortical flow is located on the flap surfaces. As a result, this higher order term can 
cause a considerable change in drag. 

On the main wing the resultant force is normal to the wing surface. From Fig. 
14 the lift and the drag components are 

L w = R w cos a (2.5.1a) 

D w = R w sin a (2.5.1 b) 

where R w is the resultant force on the main wing. 

On the flap surface there are all three components of force, as indicated by Eq. 
(2.4.7) 


R F X 

= —Rp sin e sin 8 

(2.5.2 a) 

r f v 

= Rp cos 8 

(2.5.2 b) 

r f. 

= Rp cos e sin 8 

(2.5.2c) 


where Rp is the resultant force on the flap surface as shown in Fig. 14. 
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Since the Rp z component is canceled by symmetry, only Rp x and Rp y contribute 
to lift and drag as illustrated in Fig. 14. 

Lp = Rp y cos a — Rp x sin a (2.5.3a) 

Dp = Rp y sin a + Rp x cos a (2.5.3 6) 

From Eqs. (2.5.2a) and (2.5.2b) 

Rp x = —sin e tan 8 Rp y (2.5.4) 

From Eqs. (2. 5. 2), (2. 5. 3) and (2.5.4), the drag and lift on the flap is 


Lp = Rp y (cos a + sin e sin a tan 8) 


Dp = Rp sin a (1 — 


sin e 
tan a 


tan 6 


The total lift and drag are 


L = R w cos a + Rp y (cos a + sin e sin a tan<$) 


D = sin a [ R w + Rf v 1 - 


sin e 
tan a 


tan 6 


(2.5.5a) 

(2.5.56) 

(2.5.6) 

(2.5.7) 


For small a and e these equations reduce to 

L ~ Rw + Rp y ( 1 + e a tan<5) (2.5.8) 

D ~ a (R w + Rp y (l - ^ tan 6 )) (2.5.9) 

Eqs. (2.5.8) and (2.5.9) show the influence of the forward component of force on lift 
and drag. Due to this thrust component the lift increases by an amount of order 
a 2 , while the drag is reduced by almost the same order as the drag itself. For 
moderate deflection angle, the lift-to-drag ratio is 


L_ 

D 


1 _ 

a 


1 + 




tan 8 ) + 0(a) 


(2.5.10) 


Rw + Rf v £* 

The second term within the parenthesis shows the influence of flap deflection on the 
lift-to-drag ratio. This amount is of order 1/a, the same as that of the lift-to-drag 
ratio without flap deflection. 


Chapter 3 
Straight-Feeding-Sheet Model 


In this chapter the simplest vortex model for the delta wing vortical flow, first 
suggested by Legendre and implemented by Brown and Michael, ^’ 21 is chosen to 
investigate the qualitative trends of the effect of a leading edge flap on the aerody- 
namic characteristics of a low aspect ratio delta wing. This model simply assumes 
that the rolled-up shear layer consists of a straight vortex sheet and a singular vor- 
tex core and requires this singularity system to be force free in a global sense. The 
analysis is carried out in the cross-flow plane based on the steady-unsteady anal- 
ogy developed in Section 2.2.2. Complex representation of the cross-flow plane 
is utilized. The trace of the wing in the cross-flow plane is transformed into the 
horizontal axis of a transformed plane by the Schwarz-Christoffel transformation. 
The advantage of the conformal mapping is that the boundary condition on the 
wing surface is satisfied exactly. 


3.1 Force Free Conditions 

As explained in Section 2.3.2 the force free condition is imposed on each vor- 
tex+sheet system. Combined with Eqs. (2.3.1) and (2.3.3), the force free equation 
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given by Eq. (2.3.4) for small e reduces to 

( 2a 'vj ~ °'sj) = v vj ( 3 - 1 - 1 ) 

where ()' indicates the length scale normalized with a(x), and v v j can be expressed 
as 

= £ ( F - ^1°S ('-'»;)) . (3 ' 1 ' 2) 

\ ' <7 *(J vj 

where F is the complex potential in the cross-flow plane, to be derived in the next 
section, and ( ) indicates the complex conjugate. 


In the present case there is a system of four vortices fed by straight vortex sheets 
emanating from the leading edges and from the flap hinges. Since only the sym- 
metrical configurations are treated, the analysis considers one half of the cross-flow 
plane, as shown in Fig. 15. Thus cr' s = i (k + (1 - k)e t6 ) for the leading edge 
vortex and a' s = i k for the hinge vortex. The force free conditions for both vortex 
systems are 


£(2 Ki + i(*+(i -*)«-“)) = 




(3.1.3) 


0 — y 0\)l 




+ ik ) = 


d_ 

da 



2tt 


log(<7 - a vh ) 


CT 


(3.1.4) 


where subscripts l and h refer to the leading edge 


and hinge vortices respectively. 


3.2 Schwarz-Christoffel Transformation 

The derivation of the complex potential F is facilitated by transforming the trace 
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of the wing in the cross-flow plane, together with its line of symmetry, into the real 
axis of a transformed plane. To perform this transformation the Schwarz-Christoffel 
transformation is used. The mapping procedure is shown in Fig. 16. 

The general form of Schwarz-Christoffel transformation is 

M n<0 - 0;) A (3.2.1) 

1=1 


where 0 = £ + id is the complex representation of the transformed plane. 
In integral form 


cr = A 2 + A 1 


rQ ” 

L n< e - 

•' u ° t=i 


©;/•' dO 


(3.2.2) 


where /?,- = k,-/ 7r — 1 , and K{ is the external angle of the vertex of the polygon, 
and 0 is the corresponding point of a in the mapping plane. A\ and are 
complex constants. 

From Fig. 16 


then 


7T, 

« 2 

= tt/2, 

«3 = tt-6, 

27 T, 

«5 

II 

+ 

K 6 = "w/2 


(3.2.3) 


01 =0, 02 = -1/2, = - 6 / 7T, 

04 = 1, 05 = <5/tt, 06 = -1/2 


(3.2.4) 


Riemann’s mapping theorem® 8 establishes that the correspondence of two arbi- 
trary points plus a far field condition can be imposed. Here two points are chosen 
arbitrarily as 


©4 = 0 
©o = ©2 


(3.2.5a) 

(3.2.56) 
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These conditions imply that the leading edge is transformed to the origin and the 
initial point of the integration is taken arbitrarily at 02- Requiring that the 
velocity in both planes be equal at infinity gives another condition as 

dc 


ae 


= i 


<7 — KX> 


(3.2.6) 


With Eqs. (3.2.5) and (3.2.6), the transformation function becomes 

r©' ft' /Q' - Ql\ 6 /* 


r i = r ©' f e'-e& y 

M y J (q ' - ey (©' - ey \ & - Q U 


dQ' 


(3.2.7) 


where ()' indicates the length scale normalized with a(x). 

The real constants 0^, 0g, 0g, 0g can be obtained by matching the four correspond- 
ing points in the physical and the transformed plane 


1 — k + 


1 — k 


f 


r&e 

JQ' 


• <s< 

II 

b 

-> 

©3 

(3.2.8a) 

= ik ^1 + (1 — k) e*^ 

- 

0 

(3.2.8 b) 

cr's — ik 

- 

©5 

(3.2.8c) 

4 = o 

- 

©6 

(3.2.8d) 

! a,b,c,d) to the transformation function (3.2.7) gives 

3 /0k — C 1 \^/ 7r i 

? f 5 ^ _ n 9 

\J<,% - f')K' - © 2 ) 

V0 ' -e' 

) > 

/ 

( 


6/tt 

d(' = 0 

(3.2.96) 

^(©6 - €')(€' - © 2 ) 


’* a 

/ ©5 ~ € 

\ 6/n 

) dt' = 0 

(3.2.9c) 

- (')((' ~ © 2 ) 

- 0' 


(? ~ ©5 

) de' = 0 

(3.2.9d) 

- m 1 - © 2 ) 

H' - ©3 
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This system is solved numerically for 0^, © 3 , © 5 , 0g using the Wolfe method,^ as 
described in Appendix 1. 


3.3 Mathematical Procedure 


From Fig. 16 the complex potential in the transformed plane is the sum of the 
potentials of a uniform velocity and that of four vortices, two of which are the 
image vortices needed to satisfy the boundary condition on the wing surface. Then 
the complex potential in the transformed plane is 


**-'*=» * 


and the complex velocity in the physical plane is 

dF dFdQ 


(3.3.1) 


V = 


da 


dQ da 


iTi ( 1 1 \ 

Va+ 2 *{e-e vl 


(3.3.2) 


dQ 


2tt \e-Q vh Q-Q vh J\da 


where 

d& V(9-e 6 )( e-©2) / e-e3 y /l ' 

da e \6-e 5 J *■ ' ' ; 

According to Eq. (3.3.3) the velocity becomes singular at the leading edge, 0 = 0, 
and at the flap hinge, 0 = © 5 . To avoid this singularity, a Kutta condition is 
imposed at these two points. The two equations for the Kutta condition are 

Ua + ^ (= = 0 (3.3.4) 

2tt \Q vl Q v i) 2?r \Q vh Q vh J 


Ua + 
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2tt \Q^-Q v i Q 5 -Q vl J ' 2 ?r \Q$-Q v h ©5 -&vh 




hz) - 0 (3 - 3 - 5) 
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With the complex potential F in Eq. (3.3.1), the force free conditions given by 
Eqs. (3.1.3) and (3.1.4) become 


^(2 Kl + *(* + (l-*)e-*)) = 



£l_ 1 ... + il±( — i L- Yu 

27T 27T ©i>/l ©v/ ” ®v/»/ / 
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r0v/-©3V /,r \ 

l ©,/ 

\e v i-eJ / 


+ ( iog ( 0 _ Qyi ^ ~ iog ^ ~ 


O’ *O v l 


(3.3.6) 


for the leading edge vortex system, and 



+ ik ) = 


Va + j B 


©1 >h ©1 il ®vh ~~ ©u/ 


» T /i 1 _ \ 

2 tt e vh -Q vh ) 
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V ©v/i 

\®vh~Q J / 


(3.3.7) 


+ ^“^( lo s(0-0v/>)-iog( a - <T v/ l )) o ._ 


¥ &vh 


for the hinge vortex system. 

Combining Eqs. (3.3.6) and (3.3.7) with Eqs. (3.3.4) and (3.3.5), and normalizing 
all the length scales with a(x), the following system of two complex equations for 
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the vortex positions Q v i and Q v h (or a v { and a v h ) is obtained 


e 

ka 


( 2a , vl + i (k + (l-k)e ~ i *)) = 
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(3.3.8 b) 

(3.3.9 a) 
(3.3.9 b) 
(3.3.9c) 
(3.3.9 d) 


System (3.3.8) is solved using the Newton-Raphson method, as described in Ap- 
pendix 2. 
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In the absence of flaps the only governing parameter would be the ratio of e/a. 
In the present case two additional parameters, 8 and k, are included in the group 
e/(k a ), as indicated in Eqs. (3.3.8a,b). The flap deflection angle 8 does not 
appear explicitly but through the transformation parameters ©2, ©3, ©5, ©6, which 
are functions of 8 and k . 

With the vortex positions known, the vortex strengths can be obtained from 


*T '/ 

D-B 

(3.3.10 a) 

27 r 

CB-AD 

iT' 

A-C 

(3.3.10 b) 

n 
27 r 

CB - AD a 


where T' = T/(Ua(x)) 

3.4 Aerodynamic Forces 


3.4.1 Resultant Force 

To evaluate the lift-to-drag ratio for a flapped delta wing, as derived in Section 
2.5, the pressure must be calculated first. For small a and using Eq. (2.3.1), the 
unsteady pressure coefficient given by Eq. (2.2.19b) becomes 




(?) 2 - ( t ) 2 


(3.4.1) 



3. Straight-Feeding-Sheet Model 


38 


where the v y and v z components are given by 


v y ,v 2 _ 1 dF 

W - l U - ulfo 
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(3.4.2) 


Thus 



(3.4.3) 

(3.4.4) 


where 3? and 3 s represent the real and imaginary part of a complex function. 


The v x component is best obtained from the conical expression of the velocity 
potential in Eq. (2.2.2). In terms of rectangular coordinates the potential can be 
written as 

$ c = x (U + Uay + <p c (y,z)) (3.4.5) 

where 

y = - z = - (3.4.6) 

x x 

Then the perturbation velocity in axial direction is 
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Introducing the complex potential F this becomes 


v x 

V 



O0 +— log 
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i v y 
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(3.4.8) 


The presence of the feeding sheet induces a discontinuity in the v x component of 
the velocity, leading to the treatment of branch cuts, as shown in Fig. 17. 

The resultant force is given by 

R = l p da (3.4.9) 

J each surface 


3.4.2 Normal Force 

As mentioned in the Introduction, experiments®’ 1 ®’ 1 ® have shown that the lift de- 
creases with flap deflection due to the suppression of the vortical flow. This phe- 
nomenon can be best observed by considering the lift contribution of the vortex 
systems. For small a and e the lift is related to the total force component normal 
to the wing surface 

L ~ N + Rp y ta tan <5 (3.4.10) 

where N is the total force acting normal to the free-stream and Rp y is the y- 
component of the force acting on the flap, as illustrated in Fig. 14. 

The second term represent the contribution to the total lift of the thrust component 
of force, which is of order a 2 . Thus, 


L ~ N 


(3.4.11) 
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The total normal force can be calculated by considering the momentum flux across 
the cross-flow plane 


N = - P U 


Applying Stokes’ theorem 


N 


//(I - Ua ) ** 

(3.4.12) 

= —pU y $ dz 

(3.4.13) 


The integration path is taken around the wing surface including all the singularities. 
In terms of the complex potential and the real stream function 


N = -pU $ 



i 'I' )da 


(3.4.14) 


Since the contour is a streamline 


N = -pU Of 



Substituting for F 



(3.4.15) 


(3.4.16) 


To perform the integral of the parts containing logarithmic terms, the path can be 
altered by taking a circle of infinite radius, as shown in Fig. 17, since there are no 
singularities in the far field, where dcr/dQ = 1 
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+ 2 7 log 


\0 - Q v h ) ) 




+ 
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In the transformed plane, the normal force arising from the logarithmic terms does 
not depend explicitly on the wing trace geometry, but implicitly through the vor- 
tices’ positions and strengths. This equation shows that the normal force is just 
the sum of the rate change of impulse of two pairs of vortices. Then the result of 
the integration of the logarithmic terms could be simply obtained by considering 
the impulse of a pair of vortices, or by integrating along the branch lines as shown 
in Fig. 17. Following the second approach gives 


<f P- (log (0 - ©„,) - log (0 - ©„,)) d@ 

|©| — MX) ^ 7r 


' | 0|— >00 

•r, 

2tt 


log (0 - Q v i) | p - log (0 - O vl ) | p -i 2t r dQ 


(3.4.18) 


+ / log(0 - ©„/) | p - log (0 - 0„/) 
__ JevI 

= T[ ( Qyl — 0t) /) 


dQ 


Then the normal force originating from the logarithmic terms is 

N v = pU s(r, (0„, + 0„,) + r t (0„ k -^„ft)) (3.4.19) 


This force is associated with the separated flow field about the wing. 

The component of normal force from the the first term in Eq. (3.4.16) has a 
singularity at infinity. This component could be identified with the force produced 
by the wing if the flow were attached. Integrating along the wing trace, as shown 
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in Fig. 18, there results 


Na _ l e 2 ^ e-e 5 \ </ * m 
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(3.4.20) 


where the subscript p represents the principal value (i.e., all values are calculated 
in the same Riemann surface.) 

Then the total normal force is 


N = pUZ (T,(Q vl - e vl ) + T h (Q vh - Q vh )) 

e 2 /e-Q b \ s '* 


— pU‘ 


a 5s j) 


(3.4.21) 


>/(© - © 6 )(© - ©2) 


dQ 


The normal force is composed of two parts; the normal force from the equivalent 
attached flow and that from the separated flow. 

Dividing this total normal force by the dynamic pressure pU^/2 and the total 
wing area fca(x) 2 /e, the normal force coefficient can be obtained as 

C N = \ (C Nv + C Na ) (3.4.22) 

where Cfj v and C n a represent the contribution to the normal force of the vortex 
systems and the attached-type flow respectively. Hence 

C N V ~ 


a (2r|(0;, - &,,) + 2T' h (@' vll - W vh )) (3.4.23) 
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(3.4.24) 


Here all the values of the integrand are principal values. Eq. (3.4.24) shows explic- 
itly the effect of the flap angle on the normal force through the term cos 26 in the in- 
tegration over the flap surface, while the transformation coefficients 0^, 0g, 0g, 0g 
depend implicitly on 6 and k . The region from 03 to 05 is the flap surface. 


3.5 Numerical Implementation 


3.5.1 Integration at Singular Points 


In Eqs. (3.2.9 a,b,c,d), from which the transformation coefficients 02? ©3? ©5> ©6 
are calculated, the integrands are singular at the points £ = 02, ©3, 06 • To perform 
the numerical integration, the integrand is expanded to 1st order in Taylor series 
about those points. For example, at £ = 02 


/•©2+e £ 
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202 ( ©5 - 02 

\/06 - \03 - © 2 / 


3. Straight-Feeding-Sheet Model 
where ( = £ — 02 - 


44 


Romberg integration is used on the rest of the range. ^ 

3.5.2 Integration Path for Calculating Vortex Positions 

During the numerical procedure for calculating the vortex positions in the physical 
plane through the mapping function (3.2.7), by solving Eqs. (3. 3.8a, b), taking the 
integration path as a straight line connecting ©2 with the assumed vortex position, 
as shown in the broken line in Fig. 19, causes the following problem: When the 
vortex position is very close to the real axis in the transformed plane, which is the 
case when a is small or 6 is large, the integration path passes close to the singular 
points ©3 and ©g. To avoid this possible numerical singularity, the integration 
path is taken in the following way: 

©2 — ► ©2 + i 2 — ► assumed 3? Q v + i 2 — >■ assumed ©^ (3.5.2) 

The maximum imaginary value, taken as 2, is large enough to calculate the vortex 
positions. 

3.5.3 Velocity at the Vortex Position 
Eqs. (3.3.8a,b) contain terms of the form 

^(log(0 - Q v ) - log (a - 0 -„))^_^ (3.5.3) 

which can be written as 


A( log(0 _ 0„) - log(<r - 


(3.5.4) 
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In calculating this value numerically, a finite difference scheme is employed as follows 
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(3.5.5) 


3.6 Results 

3.6.1 Conformal Mapping Parameters 

The following table shows the transformation coefficients ©[>, © 3 , © 5 , ©g for varying 
flap deflection angle with k = 0.6 and e = 15 deg. 


8 (deg.) 

©2 

©3 

©5 

©6 

40 

-0.6369 

-0.5160 

0.9785 

1.301 

32 

-0.7091 

-0.5722 

0.9531 

1.251 

24 

-0.7824 

-0.6284 

0.9183 

1.195 

16 

-0.8558 

-0.6836 

0.8794 

1.134 

8 

-0.9286 

-0.7371 

0.8362 

1.069 

2 

-0.9822 

-0.7756 

0.8010 

1.018 
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The results shows that as S — ► 0, 02 — > — 1 and ©g — ► 1, which correspond to the 
Brown and Michael case. 

3.6.2 Aerodynamic Forces; a Comparison with Experimental Results 

To compare with the experimental results of Ref. 14 a wing with t — 15 deg. and 
k = 0.6 was chosen. The angle of attack a was varied from 0 deg. to 20 deg. 

Fig. 20 and 21 show the increment of lift on the main wing and flap respectively. 
The abrupt departure from experiments observed at about 12 deg. appears to 
suggest that the trailing edge effects begin to be significant for the experiments. 

Fig. 22 depicts the total lift change. Again, the sudden departure from the 
experimental trends may indicate that those measurements were affected by the 
trailing edge effect. 

Fig. 23 shows the drag on the main wing. It can be seen that the vortex model 
overestimates the drag, particularly at higher angles of attack. 

Fig. 24 shows the drag on the flap. Notice that this component of drag can 
reach negative values. This is due to the forward force component discussed above. 

Fig. 25 shows the evolution of the resultant drag. The negative values of drag 
arise from the conical assumptions, where no trailing vortices and associated drag 
are considered. 

Fig. 26 shows the change in the L/D. No experimental values were available 
for a comparable geometry of the wing. The no-deflection case is included for 
comparison. 

Fig. 27 shows the positions of vortices for varying flap deflection. For flap 
deflection less than 18 deg. the hinge vortex almost disappears. The vortex 
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positions in the physical plane are shown in Fig. 28. It can be observed that 
that flap deflection causes the hinge vortex to move away from the surface and the 
leading edge vortex to approach the surface for the given conditions. 

Fig. 29 describes the strength of both vortices as functions of flap deflection. 
Flap deflection strengthens the hinge vortex and weakens the leading edge vortex. 
The strength of leading edge vortex reduces by a half for a 30 deg. flap deflection. 

Fig. 30 shows the different components of normal force. It is noticed that flap 
deflection acts primarily on the leading edge vortex component of the normal force. 
The hinge vortex has a negligible effect on the total force. 

Fig. 31 and 32 show the pressure distribution obtained from the straight-feeding- 
sheet model, compared with experiment.'* Experimental data are chosen at 37% 
of the chord from the apex to minimized the effects of trailing edge. Without flap 
deflection the prediction is reasonable but shows overprediction of pressure and the 
peak is shifted toward the leading edges. The pressure distribution with 15 deg. 
flap deflection angle shows unrealistic pressure peak at the leading edges. The 
unrealistic pressure peak is due to the vortex position, which is located very close 
to the wing surface. 

3.7 Summary 

Even with the crudeness of the vortex model, the two major phenomena due to 
flap deflection, increment of lift-to-drag ratio and suppression of vortical flow, are 
qualitatively observed. The leading edge vortex pair is much stronger and is re- 
sponsible for most of the lift and drag changes that take place during flap deflection. 
In agreement with experimental evidence, it is found that flap deflection causes a 
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decrease of lift and drag, and an increase of the lift-to-drag ratio. The rapid re- 
duction in drag is attributable to a thrust force component developed by the flaps, 
while the main reason for the lift reduction is the partial suppression of the leading 
edge vortical system. There is agreement with experiments in the general trends. 

This model does not yield satisfactory vortex positions, which result in inac- 
curate pressure distribution. The simplified vortex model used here leads to an 
overprediction of lift and drag forces. For better prediction of pressure distribution 
and force, a improved vortex model is required. This is the objective of the vortex 
cloud model, discussed in the next chapter. 


Chapter 4 
Vortex Cloud Model 


Due to the crudeness of the vortex model only global quantities, such as lift and 
drag, were evaluated in Chapter 3. In this chapter a vortex cloud method is adopted 
to investigate the pressure distribution and the shear layer structures. Combined 
with the steady-unsteady analogy developed in Section 2.2.2, the vortex cloud 
method leads to a numerically efficient methodology with good definition of shear 
layers. The shear layers are replaced by a cloud of discrete vortices having finite 
core size, and the core region of the vortical flow is modeled with a single vortex. 
A boundary element method is employed to represent the trace of the wing in the 
cross-flow plane. The various aspects of numerical implementation, such as core 
modeling, merging of vortices and integration methods, are investigated. 


4.1 Vortex Model 

The vortex cloud method simulating 2-D separated flow fields consists in tracing the 
trajectories of a large number of discrete vortices. A problem facing methods of this 
type is the large velocities induced by neighboring vortices. Some researchers em- 
ploy vortices with finite core size to avoid large mutually induced velocities, 60-62-65 
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while others redistribute the vortices at each time step so as to avoid logarithmically 
infinite velocities.^ In this study vortices having finite core size are employed for 
two reasons. The redistribution of the vortices at each time step is computationally 
expensive, and causes the disappearance of the coherent structures shown in Fig. 
6. This disappearance occurs because the redistribution technique prevents the 
overriding of vortices. 

The rolled-up shear layers are replaced by a cloud of discrete vortices, each of 
which is generated from the separation point at a time step. Chorin 62 first intro- 
duced the concept of vortices of finite core size in the cloud approach, so that all of 
the rotational and viscous effects can be assumed to take place in the core, while 
outside of this core the flow field is governed by potential theory. By a proper 
distribution of vorticity within the core, the infinite velocity induced by neighbor- 
ing vortices can be avoided. Using this concept, the vorticity distribution, which 
appears in the right hand side of Eq. (2.3.7), can be represented as: 

N 

«(*»*) = £r,»(l*-*i(t)l) ( 4 - L1 ) 

1=1 

where Xi is the vorticity distribution function satisfying the normalizing condition 

X»(l x l) dx = 1 (4.1.2) 

In the case of point vortices, the functions \i become the dirac-delta function 
6(|x — Xj|). Introducing a core radius q, Xi can be written as: 

X.dx-Xil) = g(!2LZB[) (4.1.3) 

where G is a shape function common to all vortices. 
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The velocity induced by the vorticity distribution r,Xj(|x — x,-(t)|) is 


u(x,<) 


27 r 


(x Xj) X ear r t - 
Ix-Xil 2 


(4.1.4) 


where 



The following vorticity distribution G and the associated velocity distribution used 
in this study were derived by Spalart^ 


7T G = 
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3(<? 2 - l) 2 


e > i 
e < i 


(4.1.6) 


2tt u <jj _ f i Q > 1 

Ti \ Q (3 - 3 e 2 + e 4 ) e < l 


(4.1.7) 


where g = |x - Xj| /q. 

It can be assumed that the most desirable core vorticity distribution would be of 
the form 1 — exp (— g ^) , since this corresponds to an isolated vortex in a viscous 
fluid. However, the computational time and effort would be considerably larger 
than in the case of distributions (4.1.6) and (4.1.7), and as shown by A. Leonard,*^ 
for fixed core size, differences in the velocity representation have a small influence 
on global flow features. The distribution used here, shown in Fig. 33, resemble 
a Gaussian distribution and was chosen as a compromise between computational 
efficiency and truthful representation of vorticity diffusion. Accordingly, Prandtl’s 
relationship 78 for core size, 

sf oc t (4.1.8) 


valid for a Gaussian distribution of vorticity, is used together with Eqs. (4.1.6) and 
(4.1.7). 
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4.2 Kutta Condition 

In this approach, a new discrete vortex is generated from the separation point 
at each time step such that the flow field including this new vortex satisfies an 
appropriate Kutta condition at the separation point. For an unsteady flow, as 
opposed to a steady flow, the vorticity is continuously supplied to the surrounding 
fluid from a separation point such that the strength and the space distribution of the 
shed vorticity allow for pressure continuity to be satisfied. 56-59 strength and 

space distribution axe functions of the velocity on each side of the vortex sheet . ^ 
From Fig. 34a, Bernoulli’s equation on each side of the vortex sheet is 

1 o dd>\ 

H\ = - p ui + pi + p — (4.2.1a) 

1 o dd> o 

#2 = - p u 2 2 + P2 + P -r£ (4.2.16) 

Here, H , u and p represent the total head, velocity and pressure respectively. The 
total head must be the same on each side since the whole flow is started from a 
uniform state. This requirement, together with pressure continuity, gives 

& -% M \ = J |(“1 2 - “2 2 )| (4.2.2) 

Since the difference in velocity potential across the sheet represents the vortex 
intensity of the sheet, there results 

y = !h 2 -“’ 2 >l 

The separation angle of the vortex sheet, 0 S , as shown in Fig. 34a, is such that the 

boundary condition on the wing surface, given by Eq. (2.3.11a), is satisfied. In 

this study the separation angle is set to zero, as shown in Fig. 34b, implying that 
vorticity going into the vortex structure is shed parallel to the lower surface. 
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This assumption makes the velocity on the upper side of the vortex sheet equal to 
zero, since the assumed streamline, as indicated in Fig. 34b, causes a stagnation 
point on the upper part of the vortex sheet. Combining Eq. (4.2.3) with the 
relationships x = Ut, the strength of a new vortex can be calculated as 


r new — 


— u\ ^ At 

2 


2 “‘ 


Ax 

IT 


(4.2.4) 


Normalizing with Ua(x ), Eq. (4.2.4) reduces to 

w __ F new 1^ /u \ \ 2 Ax k 

new Ua(x ) 2 V U / x tane 


(4.2.5) 


The position of this new vortex can be obtained by considering the convection 
velocity at the separation point. By definition, a vortex sheet moves with the 
average of the velocities at its two sides. Then, the convection velocity of the 
newly shed vortex is 

u s = ^( ui+U 2 ) = ^ U 1 (4.2.6) 

The position of the new vortex is then 

A x 

<7new = @8 ~t" 6$ At Ug — <7 g -f- Ug (4.2.7) 


where a s is the position of the separation point, and e s is the unit vector tangent 
to the lower surface at the separation point, as illustrated in Fig. 34c. 

Normalizing with a(x ) , Eq. (4.2.7) becomes 
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(4.2.8) 
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4.3 Wing Cross-Section Representation 

The problem consists in solving Eq. (2.3.7) with the vorticity distribution given by 
Eq. (4.1.1) on the rolled-up shear layer, subject to boundary conditions (2. 3. 11a, b) 
and the Kutta condition expressed by Eqs. (4.2.5) and (4.2.8). 

As mentioned in Section 2.1, the boundary value problem given by Eqs. (2.3.7) 
and (2.3.11) can be solved in terms of a distribution of sources and doublets on the 
boundaries. The doublet distribution is replaceable by a vorticity distribution.®* 
In this study the wing cross-section is represented by vorticity panels of linearly 
varying strength. 

If the body is divided into m — 1 vorticity panels and there are n vortices in the 
shear layers, as shown in Fig. 35, the stream function at the point i in the flow 
field can be written as 


u(aj) K(oi, Gj) da j 

1 " 

+ ^~J2 r k lo S r tJfc + U tan a yi 


(4.3.1) 


k=l 


where w(aj) is the bound vorticity density at the j ^ panel, and K(a{,aj ) are the 
geometric influence coefficients, which depends only on wing geometry. r,j. is the 
distance between the i point and the vortex. 

The first term on the right-hand side is the stream function due to the bound 
vorticity panels, and the second term represents the stream function due to the 
discrete vortices distributed on the shear layer. The last term is due to the vertical 
component of the free-stream. 
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Boundary condition (2.3.11a) gives 


$1 = #2 = 


= 'I'm = const. = c 


(4.3.2) 


Conservation of the total circulation in the flow field as given by Eq. (2.3.9) 
implies that the total vorticity in the whole flow field must be zero 


/ 


ujjdaj + 5Z r ib = 0 
Jfc=l 


(4.3.3) 


Here there are m + 1 unknowns; m u's for the vorticity panels and c, the value 
of the stream function on the wing trace in the cross-flow plane. Eqs. (4.3.1) and 
(4.3.3) are a set of m + 1 equations. Combining Eqs. (4.3.1), (4.3.2) and (4.3.3), 
the following system of equations results 
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The matrix a^j is called the geometric influence coefficient matrix and is constant 
during time iteration. The bj ’s are the stream function values induced by the vor- 
tices in the shear layers and the free stream. The last row in the matrix represents 
the integration of the bound vorticity panels. Using the trapezoidal rule for the 
integration, dj can be written as 

dj = \dj + lj+ 1) (4.3.5) 

where lj is the length of vorticity panel. 
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4.4 Solution Procedure 

The steady-unsteady analogy developed in Section 2.2.2 shows that the solution of 
Eq. (4.3.4) (bound vortex strength and stream function at wing trace) as well as the 
the position and strength of the shear layer core (represented by a single vortex) 
must reach a converged state when expressed in terms of a similarity variable. 
Converting the time variable to the x coordinate through the relationship x = Ut, 
all the length scales can be normalized with a(x) and the vortex strength with 
U a(x ) . Since in this transformation the total wing span length a(x) increases with 
time, the vortex strength, vortex position and core radius are re-normalized at the 
beginning of each time step. Only the right-hand side vector of Eq. (4.3.4) is varied 
at each time step. Such changes result from an increment in the number of vortices 
in the shear layer through the Kutta condition given by Eqs. (4.2.5) and (4.2.8), 
and from the movement of the vortices as a consequence of the induced velocities. 
The left-hand side matrix remains constant throughout the calculation procedure. 

4.5 Aerodynamic Forces 

The pressure coefficient was derived in the previous chapter as 

C P = tan’a - 2^ - (£)* - (£) 2 (4.5.1) 

The magnitude of the velocity in the cross-flow plane is numerically equal to the 
bound vortex strength per unit length, since the flow internal to the wing cross- 
section is assumed to be at rest.** 7 As indicated in Fig. 36, the bound vortex 
strength can be written as 
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hence 

oj = ^ Vy 2 + v z 2 (4.5.26) 

In this study, the bound vortex strength itself represents the outer velocity of the 
boundary layer. As a result, the magnitude of the velocity in the cross-flow plane 
can be obtained from the solutions of Eq. (4.3.4). The perturbation velocity in the 
axial direction can be obtained from the following conical relationship 


v x — — (U tana + (f> u — y v y — z v z ) 


(4.5.3) 


The velocity potential is given in terms of logarithmic functions and its evaluation 
requires consideration of branch cuts, as shown in Fig. 37. By integrating the load 
distribution given by Eq. (4.5.1) over wing and flap surfaces, the total resultant 
force on each surface can be obtained. 


4.6 Trailing Edge Wake 

The trailing edge wake of a delta wing, as shown in Fig. 1, has a shape different 
from that of a conventional large aspect-ratio wing. The pressure peak just below 
the leading edge vortex core reflects a vorticity distribution such that the bound 
vortex strength, as a function of spanwise location, has a maximum between the 
wing center and the wing leading edge, as illustrated in Fig. 38. Conical flow 
implies the assumption that the wing extends infinitely far downstream, namely, 
that no trailing edge exists. To study the trailing edge wake roll-up in a qualitative 
way, an artificial trailing edge is assumed to exist, from which the conical vorticity 
distribution calculated on the wing surface is released. A Kutta condition is applied 
such that there is no pressure jump at the trailing edge. 

The bound vorticity on the wing surface has two components, as shown in Fig. 
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39. These two components must satisfy the divergence free condition 


_ =* du> 

V • Q. — -jr— + 
ox 


dr] 

Tz 


0 


(4.6.1) 


where Q = u e x + V e z 

The absence of a pressure jump at the wing trailing edge demands that the shed vor- 
ticity at the trailing edge should be aligned with the local velocity. An approximate 
form of the trailing edge Kutta condition 92 can be written as follows 


d rj 

7Tz 


0 


then Eq. (4.6.1) reduces to 



(4.6.2) 


(4.6.3a) 


or 

oj = u>(z) (4.6.36) 

This indicates that uj, which is identified with the w’s from Eq. (4.3.4), is constant 
as it moves downstream. After the trailing edge, the evolution of the wake is solved 
in a straightforward way using a 2-D time evolution technique. 

4.7 Numerical Implementation 
4.7.1 Influence Coefficients 

The stream function at a control point i due to the panel, having linearly 
varying vorticity intensity, and its corresponding panel at the other side of the 
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plane of symmetry, as shown in Fig. 40, can be written as 


(^j + J ~ ) (logrjs — logrjj) ds 


= p ijUj + Qij + 1 (wj+i - uj) 

where lj is the length of j th panel and 
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Evaluating these integrals, Pij and Qij can be written as 
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The stream function value at point i due to the bound panels is 
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(4.7.1) 

(4.7.2a) 

(4.7.26) 

(4.7.2c) 

(4.7. 2d) 

(4.7.3a) 

(4.7.36) 


(4.7.4) 
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Then the elements of the matrix of influence coefficients are 


a 


ij 


— 2jr (^*j Qij + Qij- 1) 


(4.7.5) 


4.7.2 Induced Velocity 


Once strength of the bound vortex panels is obtained by solving Eq. (4.3.4), the 
flow field velocity is used to update the positions of the vortices in the shear layers. 
From Eq. (4.1.4), each component of the induced velocity outside the core is 


v v = 
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Z — Zj 
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Vz = -7T- 


2?r (V-Vj) + { z - z jY 
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(4.7.6a) 

(4.7.66) 


2?r (V-Vj) + (z-ZjY 
where Tj is the strength of the j th vortex. 

With 

V ~ Vs = V - {yj + syj) (4-7.7 a) 

z — z s = z — (zj + szj ) (4.7.76) 

where ( ilj,Zj ) is the slope of the vortex panel, the velocity induced by the 
bound vortex panel can be obtained from Eqs. (4.7.6a, b) with the following 
definition of T j 

Tj = (u >j + 5 ) ds (4.7.8) 

When the point where the induced velocity to be calculated is located within the 
core, the induced velocity is obtained from the velocity distribution in Eq. (4.1.7) 

(4.7.9a) 


r,- , , 3 - 3 g 1 + e 4 


V V = d (z - Zi) 


.2 


Tj 3 - 3 q 2 + e 4 

= - 27 W ~ 


.2 


(4.7.96) 
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where 


e 2 = 


(y ~ Vj? + (z - zjY 




x 


(4.7.10) 


4.7.3 Numerical Implementation of the Kutta Condition 

Flow separation in a 2-D or a 3-D flow arises from viscous effects. An equivalent 
effect is imposed on the inviscid problem through a Kutta condition. A converged 
quasi-steady analysis requires that the shed vorticity at the wing side edges be 
convected with finite velocity. This velocity, not known a priori, is the average of 
the velocity on both sides of the wing at the side edge. ^ 9 In a 2-D case the velocity 
at the separation point can be determined from the experimentally established 
fact that the Strouhal number, when expressed in terms of the base pressure at 
the separation point, is constant. 93 Experimental observations on the flow past an 
inclined flat plate have shown that this velocity is about 1.5 U - 94 Under the conical 
flow assumption, the value of this velocity must be constant along the leading edges. 
Other workers have attempted to find this value experimentally and have applied it 
to the solution of delta wing problems with the vortex cloud method, 31 failing to 
attain converged solutions. In this work this velocity is theoretically determined 
by recognizing the fact that, after a sufficiently large time and for given values of 
a , e and 6 , it must achieve a constant value. It was found that if an initial value 
for the bound vortex strength of about 317 is assumed, convergence to the final 
value of the side edge convective velocity is achieved within 300 steps. Too large 
initial values lead to oscillatory end behavior, while too small initial values require 
significantly more time steps. 
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4.7.4 Core Model 

The core region for a conical flow has large axial velocity, and the flux of vorticity 
generated along the leading edges and supplied to the core through the rolled-up 
shear layer is balanced by that convected by the axial velocity. Thus, the core 
region, usually identified with the “vortex jungle” in a 2-D vortex cloud computa- 
tion, cannot be included in this model of the conical flow since it would violate the 
condition of aerodynamic slenderness. Hence, the core region is represented by a 
single vortex through a simple merging scheme. By imposing a restriction on the 
rotation angle of the vortices, the efficiency of the numerical effort is increased and 
good definition of shear layers is achieved. Other studies ^>37 h ave shown that the 
rotation angle of the shear layer, as long as it is set greater than 2n, has a small 
effect on overall quantities. The same conclusion is arrived at in this study. 

The vortices which rotate more than a given angle are allowed to merge into 
the core, which is then simulated by a single vortex. The position, strength and 
core size are determined such that the angular momentum before and after the 
merging is conserved. If k vortices are merged together, conservation of angular 
momentum outside the core region is expressed by the following relationships for 
the core position and strength 

k 

E *i r i 

a c = 2=1 (4.7.11) 

E r, 

i= 1 
k 

r c = £ r.- 

t=l 


(4.7.12) 
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The angular momentum within the core is 

M = f f (<7 c + r) u r d6 dr (4.7.13) 

JO Jo 

where a c is the position of the core center, and r is a distance from the core center 
to a point inside the core, u is the velocity distribution inside the core as given by 
Eq. (4.1.7). Performing the integration, the angular momentum reduces to 

M = (4.7.14) 

Then the core size after merging is 


Cc 




E r ,■ q 3 
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V E Ti 

\ i= 1 


(4.7.15) 


Strength and position are independent of the vorticity distribution inside the core, 
but core size is dependent on vorticity distribution. 


4.7.5 Merging between Neighboring Vortices 

A similar merging technique as the one used for the core model is also applied 
to the neighboring vortices such that if the distance between them is less than a 
given value two vortices are allowed to merge into one. This value is set equal to 
the mean value of the radii of the two vortices. Position, strength and the core 
radius after merging are determined from Eqs. (4.7.11), (4.7.12) and (4.7.15). The 
merging scheme allows for a reduction of the total number of vortices by almost a 
half, making the scheme much more efficient. 

4.7.6 Absorption of Vortices on Wing Surfaces 


In this discretization of a separated wake, vortices coming very close to the surface 
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have the tendency to penetrate the wing surface. This is observed to happen even 
with very small time step and large number of panels. A means of resolving this 
difficulty has been suggested by Lighthill,®® and consists in allowing the vortices 
that come within a thin layer on the wing surface to disappear. Another way 
consists in relocating the vortices which penetrated the surface in a preceding time 
step, to its image point outside the surface. Here the vortices which come closer 
than a given distance to the wing surface are made to disappear. This distance is 
set equal to 1% of the local wing span. This surface merging tends to occur when 
the angle of attack is small or when the flap deflection angle is large. 

4.7.7 Integration Method for Advancing Vortex Position 

Any 2-D discrete vortex method contains a vorticity diffusion effect due to the 
numerical error involved in time integration. The numerical diffusion effect is 
shown in Fig. 41. Some authors®®’®® maintain that this numerical error could 
be responsible for instability of the vortex motion. Others®®’ 7 ® suggest that this 
numerical error simulates a turbulent diffusion rather than a viscous diffusion. In 
a conical flow the diffusion effect arising from the numerical error is less marked 
than in a 2-D flow. This may be attributed to the fact that the shear layers in 
the conical flow develop in a field of intense circumferential velocity, which allows 
for a shorter time during which numerical dissipation acts. Here both a 1st order 
Euler method and a 4th order Runge-Kutta method were used for time integration. 
These schemes are: 

1st order Euler method 


<T v (t + At) = a v (t) + v v (t, <r v (t)) At 


(4.7.16) 
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4th order Runge-Kutta method 

a * (* + ' t ) = av ® + Y Vv ^’ <Tt; ( < )) 

< (* + = ^(<) + -y Vt) ( t + ( <+ ^”)) 

*;**(< + At) = a,(t) + Aiv, (f+^,<r**(t + At/2)) 

(^(t+Af) = <7„(*) + ^ V v (t,cr v (t )) + 2v v |/+y,^<+yJJ 

+ 2 v v + — , <7** ^ ^ + v„ (t + At, cr***(t + At)) 

(4.7.17) 

4.8 Results 

4.8.1 Shear Layer Geometry and Surface Pressure Distribution 

Fig. 42 shows the convergence of the vorticity intensity at the separation point 
for an initial guess of 317. Significant oscillation can be seen during the first 50 
steps. If too large an initial guess is used, the vortex strength oscillates such that 
it changes sign during some time steps, leading to nonphysical results. Given a 
and 8 , a proper initial guess for each case is needed to obtain the final vorticity 
intensity at the separation point. This value decreases as the flap deflection angle 
increases and as the angle of attack decreases. 


The convergence of the present calculation can be assessed through core position, 
core vortex strength and core radius, as shown in Figs. 43a, b. The core vortex 
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strength shows some oscillation even after convergence is attained, since the vortices 
in the shear layer merge into the core at finite intervals. 

Fig. 44a shows the vortex jungle that arises when no core model is used. In 
addition, the shear layer is not clearly defined. Fig. 44b shows a case where a core 
rotation angle of 2.57 r was imposed for the same configuration and exhibits good 
shear layer definition. 

The effect of the merging scheme between neighboring vortices is shown in Figs. 
45a, b. By replacing high density clouds with fewer vortices, the total number of 
vortices is reduced by almost a half. 

Fig. 46 shows the effect of two schemes of time integration. As anticipated in 
Section 4.7.7, the 1 st order Euler scheme gives a much more diffused shear layer 
than the 4^ order Runge-Kutta method. Overall flow field features do not depend 
significantly on the integration scheme. 

The data for comparison with experiments were taken from Ref. 3, where the 
model had a trailing edge thickness of 7.7% of the half span and a 22 deg. half 
apex angle with undeflected flap. The ratio of the main wing span to the total 
span was 0.6. The data used for comparison were measured at a spanwise station 
located 37% of the chord from the apex, to minimize the trailing edge effect. 

Figs. 47a, b show the rolled-up shear layer and surface pressure distribution for 
a =25 deg. and 6=0 deg. Agreement can be considered very good, except for a 
small region near the leading edge, where a secondary vortex may have been present 
in the experiments and would be responsible for the disagreement. 

Figs. 48 and 49 are the cases for flap deflection angle of 15 deg. and 30 deg. 
and same angle of attack as the case of Fig. 47. The pressure at the centerline 
decreases with flap angle, which indicates that the effect of flap deflection is similar 
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to a reduction of angle of attack. Although the vortex position does not move much 
with flap angle, some movement toward the surface is observed, which has also been 
reported elsewhere. 1 ^ The hinge vortex is hardly observable in the cases analyzed 
here, being quite weak and tending to be absorbed into the wing surface. Fig. 49 
shows a weak hinge vortex. The greater discrepancy between computation and 
experiment near the leading edge is attributed to the stronger secondary vortices 
at higher flap deflection. 

Fig. 50 shows the case of a = 15 deg., without flap deflection. The pressure 
distribution shows good agreement. The effect of a secondary vortex is also stronger 
than at a =25 deg. The width of the pressure peak region is less than shown by 
the measurements. 

Fig. 51 shows the case with 15 deg. flap deflection angle with a = 15 deg. This 
figure shows the same trends as the straight-feeding-model, shown in Fig. 32. 

Fig. 52 shows the rolled-up shear layer and surface pressure distribution for a — 
35 deg., without flap deflection. The shear layer shows a tendency to produce the 
type of coherent structures which have been experimentally observed, as described 
in Fig. 6, taken from Ref. 77. 

The cases of 6 =15 deg. and 35 deg. with same angle of attack are shown in 
Figs. 53 and 54. In both cases, the hinge vortices have almost entirely disappeared. 
For a flap deflection of 15 deg. the pressure peak is underpredicted. The effect of 
the secondary vortex is clearly seen in the case of 35 deg. flap angle. Calculated 
results show a strong adverse pressure gradient near the leading edge on the upper 
surface, as indicated in Fig. 54, which would cause flow separation. This separated 
flow, not resolved for in this calculation, would give rise to a strong vortical flow, 
which results in the appearance of a secondary pressure peak. This secondary peak 
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is clearly seen in the experimental data in Fig. 54. 

In Fig. 55, the calculated trailing edge wake is compared with computations from 
a higher-order panel method. Reasonably good agreement is obtained. 

4.8.2 Lift, Drag, Lift-to-Drag Ratio 

The following table shows computed and measured changes in each of the force 
components and the lift-to-drag ratio as a result of flap deflection, for a = 25 deg. 
The experimental data are sectional force coefficients taken at 37% of the chord 
from the apex. 
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While lift and drag decrease with flap deflection, the lift-to-drag ratio increases. 
For a flap deflection of 30 deg., that increment is about 30%. 


4.9 Summary 

The vortex cloud model gives satisfactory agreement with experiments in pressure 
distribution and global forces. Good shear layer definition is also obtained. 

The wing model analyzed in this work exhibited negligibly small hinge vortices, 
both in theory and experiment. The usually observed instabilities of delta wing 
shear layers are hardly found, except in the high angle of attack range, and become 
more visible if no merging scheme between neighboring vortices is implemented. 
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Use of a 1 st order Euler scheme for time integration shows more clearly this kind 
of instability than a4 t * 1 order Runge-Kutta method. 

Comparison of pressure distributions with experimental measurements shows 
that the effect of the secondary vortex becomes significant in the range of small 
angle of attack and large flap deflection angle. Limitations of this model can 
be observed through comparison with experimental results for large flap deflection 
angles. The disagreement is caused partly by the effects of secondary vortices, 
and partly by the inability of the cross-flow plane analysis to resolve true velocity 
components on the flap surface. 

Leading edge flap deflection reduces the lift and drag of the wing and increases its 
lift-to-drag ratio. As was the case with the simpler model in the previous chapter, 
it is found that flap deflection is equivalent to a reduction in angle of attack. 




Chapter 5 
Comparison of Results 


In this chapter the results obtained from the two approaches analyzed in Chapter 
3 and 4 are compared with results from other analyses and experiments. Exper- 
imental data are obtained from Ref. 3, and other analyses are given by Smith’s 
theory 30 and a 3-D panel approach. 19 

Fig. 56 shows the comparison of core position and rolled-up shear layer geometry 
between Smith’s calculation and the vortex cloud model, for a/e = 1. Smith 
determined the shear layer core position by applying the force free condition to the 
core, which was modeled as a concentrated vortex. In this respect his core model 
differs from the vortex cloud’s, which uses a simple merging scheme by imposing a 
restriction on the shear layer rotation angle. Although the core models are quite 
different, shear layer shape and core position show good agreement between both 
methods. 

A comparison between results from a 3-D panel analysis and the vortex cloud 
model is illustrated in Fig. 57 for the case of a leading edge fence, corresponding 
to an upward deflection of the leading edge flap of more than 90 deg. The results 
from the 3-D panel analysis were taken from Ref. 19 for e = 16 deg. and k = 0.83, 
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and chosen at the 30% of the chord from the apex to minimize the upstream effect 
of the trailing edge. Fig. 57a shows the pressure distribution without the fence for 
a = 20 deg. The pressure peak region calculated with the vortex cloud model lies 
closer to the leading edge than that from the panel calculation. The case of 130 
deg. fence angle is shown in Figs. 57b, c for the same angle of attack as Fig. 57a. 
The results show better agreement than without fence deflection. 

A detailed discussion of the computation between theoretical results from this 
work and experiments is carried out, and possible reasons for the disagreements 
with experiment are discussed. 

Pressure distributions for a = 25 deg. and <5 = 0 deg. are shown in Fig. 58. 
Also shown are the results from the straight -feeding-sheet model. As mentioned in 
Chapter 3, the straight-feeding-sheet model overpredicts the pressure peak, which 
appears shifted toward the leading edge, while the pressures on the lower surface and 
in the midspan region on upper surface are well predicted. This implies that the 
straight-feeding-sheet model seems to primarily overpredict the vortex component 
of lift. The results from the vortex cloud model agree well with experiments. 

The pressure distribution for the case of a = 25 deg. and <5 = 15 deg. is shown 
in Fig. 59. The results from the vortex cloud model show that the pressure on 
the upper surface is almost constant near the center region and rises more sharply 
to a pressure peak than the experimental results. A similar situation occurs in 
the region where the vortices are absorbed by the wing surface. This phenomenon 
can be more clearly seen in case of <5 = 30 deg., as illustrated in Fig. 49. The 
overall pressure distribution agrees well with experiments, while the pressure peak 
is slightly overpredicted. In this case, the hinge vortices are extremely weak. This 
fact was also suggested by results from the straight-feeding-sheet model. The 
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extreme closeness of the hinge vortex to the wing surface caused oscillations in the 
numerics of the straight-feeding-sheet model, such that it was found convenient to 
neglect the hinge vortex altogether. Results obtained neglecting the hinge vortex 
are shown in Fig. 59. The pressure peak is unrealistically overpredicted and the 
pressure peak region is much narrower than indicated by the vortex cloud model or 
experiments. This is explained by the position of the leading edge vortex being very 
close to the wing surface. 9 ® Thus, the straight-feeding-sheet model exaggerates the 
vortex displacement with flap deflection. 



Chapter 6 

Conclusions and Recommendations 


The effects of leading edge flaps on the aerodynamics of delta wings was investigated 
theoretically using two different approaches. A mathematical analogy was invoked 
to treat a 3-D conical flow field in subsonic and supersonic cases. Through the 
analogy, a conical steady flow field can be represented by a 2-D unsteady self- 
similar ideal flow, whose treatment is far simpler than that of a 3-D flow. Two 
different vortex models were used as techniques to solve a 2-D unsteady flow, and 
the results were compared with experimental data and other methods of analysis. 


6.1 Conclusions 

6.1.1 Straight-Feeding-Sheet Model 

While this simple model leads to poor agreement in pressure distribution with exper- 
iments, some of the the effects of leading edge flaps on the delta wing aerodynamics 
are suggested. The wing’s geometric parameters, such as e, 6 and &, together with 
the angle of attack appear explicitly in the equations for the vortex positions and 
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the normal force components. From Eq. (3.4.23), which shows the normal force 
arising from the vortex system, it can be predicted that the normal force component 
reduces with flap deflection, since the leading edge vortex strength and the distance 
between the vortex and its image vortex in the transformed plane reduce with flap 
deflection. Also Eq. (3.4.24), which is the normal force component corresponding to 
an attached flow situation, shows the explicit effects of flap deflection, resulting in a 
decrement of the normal force component with flap deflection. This model indicates 
that the increment of lift-to-drag ratio during flap deflection is more significant in 
the low angle of attack range, as observed in experimental results. ^ The straight- 
feeding-sheet model exaggerates the vortex displacement and leading edge suction 
arising from flap deflection. 

6.1.2 Vortex Cloud Model 

This model provides a good definition of shear layers and good agreement in surface 
pressure distribution with experiments and 3-D panel calculations, as long as the 
validity criteria derived in Chapter 2 are respected. The effects of flap deflection are 
properly represented. The non-appearance of hinge vortices in the range of angle of 
attack analyzed in this approach is confirmed by experimental observation. Cases 
of low angle of attack and large flap deflection angle would violate the conicality 
and slenderness criteria and are consequently not expected to be predicted well by 
this model. 

6.1.3 Effects of Flap Deflection on Delta Wing Aerodynamics 

A decrement in vortex strength, caused by flap deflection, delays vortex breakdown 
by reducing the adverse pressure gradient around the trailing edge. Lift and drag 
decrease with flap deflection, while the lift-to-drag ratio increases. The reduction 
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of lift and drag arises from the partial suppression of the vortical flow, and the 
increment of the lift-to-drag ratio derives from the thrust component of force acting 
on the flap surfaces. 


6.2 Recommendations 

The vortex cloud model gives good definition of shear layers, which current Navier- 
Stokes solvers cannot provide since an impractically fine grid distribution would be 
required to capture such discontinuity surfaces. This work can provide some insight 
in supplying initial grid distributions to Navier-Stokes solvers, thereby reducing 
computational cost. 

This study could also be extended to moderately supersonic flows where the Mach 
cone angle is closer to the wing apex angle. In such cases the approach as developed 
here would fail, since the far field boundary conditions would be altered by the Mach 
cone. A mapping function, the so-called Tschaplygin transformation, 82 has been 
developed to account for the Mach cone trace in the cross-flow plane. Combined 
with the Tchaplygin transformation, the steady-unsteady analogy developed in this 
study would provide a scheme to treat a moderately supersonic conical flow. 

In this study the comparison between the results from experiments and calcu- 
lation based on the steady- unsteady analogy is made only for the incompressible 
case. For compressible flows, application of a similarity rule for the pressure co- 
efficient would be required, since the wing geometry and the free-stream condition 
are distorted as a result of applying the Prandtl-Glauert transformation. Such a 
similarity rule for the pressure coefficient for a slender wing is given in Appendix 
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(a) Plane view 


leading edge vortex 



(b) View from the trailing edge 



Fig. 1 Overall flow field past a plane delta wing 
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Tangential blowing: 

Stabilizes vortical flow 
at high angle of attack 

(a) Blowing (b) Leading Edge Flap 


Fig. 2 Concepts for controlling the vortical flow 



Fig. 3 Simplified flow field past a flapped delta wing 
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Fig. 4 Comparison of L/D for three different wing configurations, taken form Ref. 2 






Fig. 6 Characteristics of the shear layer structure of leading edge vortices 
compared with the plane mixing layer 
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(a) 3-D Conical vortex sheet 



Fig. 7 Conical and 2-D unsteady vortex sheets 
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singular vortex core 



Fig. 9 Straight-feeding-sheet model 



Fig. 10 Integration contours for the definition of circulation 
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(b) Transformed plane 


Fig. 16 Mapping procedure 



Fig. 18 Integration path along the wing trace 
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Fig. 20 Lift increment on main wing, 6 = 40 deg. 
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Fig. 23 Drag on main wing, 6 = 40 deg. 
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Fig. 24 Drag on flap, 8 = 40 deg. 
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Fig. 25 Total drag, 6 = 40 deg. 





Figures 


106 














109 



Fig. 32 Surface pressure distribution for a = 25 deg., 6 = 15 deg. 
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(a) Velocity distribution 
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(b) Vorticity distribution 
Fig. 33 Spalart’s vortex model 
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Fig. 36 Boundary layer and outer velocity 
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(b) For wake vortices 


Fig. 37 Branch lines 
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Fig. 38 Chordwise vorticity distribution 
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Fig. 42 Convergence of the vorticity intensity at the separation point, initial guess is 3i7 
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(b) With core model, 2.5n rotation angle 


Fig. 44 Effects of shear layer core model, a = 25 deg. 
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(b) With merging scheme 


Fig. 45 Effects of merging scheme between neighboring vortices, a = 25 deg. 
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Fig. 46 Effects of integration scheme on shear layer geometry, a = 35 deg. 
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Fig. 53 Shear layer geometry and surface pressure distribution, 
a = 35 deg., 6 = 15 deg. 
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(a) Vortex cloud model 


(b) Higher order panel method, 
taken form Ref. 61. 



Fig. 55 Trailing edge wake evolution for a plane delta wing 
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(c) Fence surface for 6 = —130 deg. 


Fig. 57 Comparison with 3-D panel analysis 
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Fig. 59 Surface pressure distribution, a = 25 deg., 8 = 15 deg. 
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Appendix 1. 

Wolfe Method 

Eqs. (3.2.9a,b,c,d) are solved numerically using Wolfe method. Wolfe’s method is 
a variation of the secant method for a system of nonlinear equations. This method 
has the advantage that as the solution becomes closer to the root, the convergence 
speed is much faster than that of Newton- Raphson method.*^ For a single equation 
the secant method is equivalent to solving the following two equations, with two 
initial guesses 


PI + P2 = 1 
Pl/Ol) + P2/O2) = 0 


The solution is updated as 


*3 +1 = PI z? + P2 x 2 


(A. 1.1a) 

(A.l.lb) 


(A- 1-2) 


For a system of n equations a set of n + 1 equations must be solved with n + 1 
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vectors of initial guesses, 

n+1 

izpj = 1 

j~ i 

n+l 

^ y Pj ■ ‘ ■ ) X rx) — 0 1 = 1^2,- • ■ , 7T. 

j = 1 

For the solution update we ignore the vector x J , which gives the maximum residue. 

n+1 

x ^lfc+l tA step — VI Pi x ^'ljb tk step (A. 1.4) 

j= 1 

where = (zj, x^, • • • , x J n ) 

This procedure is repeated until the maximum residue becomes smaller than a given 
value. 

Appendix 2. 

Newton-Raphson Method for Complex Functions 

Eqs. (3.3.8a,b) are solved using the Newton-Raphson Method. Since these are 
complex equations, the conventional Newton-Raphson method is modified accord- 
ingly. 

For the system of real equations 

/t(x) =0 i = 1,2, •••,« (A.2.1) 

the general form of the Newton-Raphson method is 

x„+i = x„ - J~ l /f(x„) 


(A. 1.3a) 
(A. 1.36) 


(A.2.2) 
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where J is the Jacobian of /i(x) at x n . 

In the case of a system of complex equations 

x = (yi+izi,y 2 + iz 2 ,---,y n + izn) (A.2.3) 


Applying the scheme to the real and imaginary parts of f , there results 


with 


J 


( 


3 


5 


( y i\ 

*i 

V2 

z 2 


/yi\ 
Z 1 
V2 
z 2 

- J- 1 

fS f j\ 

3 fl 
3® /2 
3/2 


n+1 


n 

*/» 

\3/n/ 

$n) 

*(£) 



*(U 

ir) 

®(&) 



HU 



... 

*(b) 



*>(&) 

... 


a® 


(A. 2.4) 


(A. 2.5) 


For the present problem there are two complex equations for the vortex positions, 
resulting in a 4 x 4 matrix. The reason why the Newton-Raphson method is 
preferred over the Wolfe method in this case, is that the latter needs a very accurate 
initial guess in order to converge. 


Appendix 3. 

Similarity Rule for Pressure Coefficient in Compressible Flow 

Gothert’s similarity analysis^ is modified and applied to a general slender body at 
angle of attack. 
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The pressure coefficient for a compressible flow using the small disturbance as- 
sumption is 


c p 



(‘ 



( 30 2 - Cwf 


(A. 3.1) 


The analysis consists in considering two flow fields characterized by Mach num- 
bers M\ and M 2 , with body shapes identified by the functions F\ and . Using a 
cylindrical coordinate system, the small disturbance potential equation for the flow 
field having free-stream Mach number Mi and velocity U\ with incidence angle ot\ 
can be written 


d 2 <j> 1 1 ( d 2 (f> 1 ldfa 1 d 2 <f>\ \ _ 

dx ^ 1 — Mi^ \ dr ^ r dr r 2 dd 2 / 

The boundary shape is given by 


(A.3.2) 


F\ = r — b\c fi(x/c,6) = 0 (A.3.3) 

where b and c represent the aspect ratio and chord length of the wing respectively. 
The boundary condition on the wing surface given by Eq. (2.1.7) can be written as 


V^i|,-VFi = 0 (A.3.4) 

where s indicates that V</q is evaluated at the body surface. 

Using Eqs. (A.3.2) and (A. 3. 3), Eq. (A. 3. 4) reduces to 


- tfi + 


d<f> 1 
dx 


\bi c ^ + fui tan«i sin0 + ^- 


1 ( TT l a , d<f> 1 ^ dfi 

~ { Ult ^ aicose+ w.)w = 0 


(A. 3.5) 


Invoking the slender body assumption® 1 


d<!> d<f> 1 d<(> 

dx < dr and rd6 


(A.3.6) 
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Eq. (A. 3.5) reduces to 


d<t> l 


dr 


dh 


= — U\ tanai sin# + U\b\c 
s ox 


+ 7T ( U\ tanaj cos 9 + 

heft \ 


dh_ \ dfi 
89 s) 89 


(A. 3.7) 


and the pressure coefficient given by Eq. (A. 3.1) becomes 


C P 1 = ~ rr- 


2 8<f>i 


Ui 8x 


S U 1 


1 ( ^^1 \ ^ 

— 2 ( U\ tanai sin# + — j 


- I (Pl‘-«lcos«^)_ 


(A. 3.8) 


Next consider another flow field having Mach number M 2 and velocity U 2 with 
incidence angle 02 . The velocity potentials of the two flow fields are related as 
follows 

^ = A ^ 2 ( X, \Zl -M| r ’^) (A3 ' 9 ^ 

where A is a constant to be derived later. 

The new coordinates are 

X = x, f = y 1 _ M \ r, 9 = 9 (A.3.10) 

Substituting Eqs. (A. 3. 9) and (A.3.10) into Eq. (A. 3. 2) the latter reduces to 


d 2 <fc , 1 ( d 2 (j>2 1 dfo 1 d 2 <f>2 \ _ n 

8x 2 1 — M 2 2 V dr 2 f dr f 2 89 2 / 

The body shape for the second flow field is defined by 


(A.3.11) 


F 2 = r - &2 C f2(x/c,9) 


(A.3.12) 
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Accordingly, the boundary condition is 

d<t>2 


dr 


df 2 


= — U 2 tan «2 sin $ + U 2 b 2 c — 

ox 


hcfi 

and the pressure coefficient 

n 2 d<f>2 

C P 2 = - 77 - 


+ 7—737 ( U 2 tan 02 cos# + 
cfo \ 


d<t>2 

ae 


\dh 
*) d~e 


(A.3.13) 


u , 2 ( 


~ dd>o \ ^ 

U 2 tan 02 sin 9 + -q~-J 


U 2 dx 

772 7 TT75 ( ^2 tan a 2 cos 6 + 

£V h. l c l $ 2 


d<f>2 V 
ae ), 


(A. 3.14) 


Using the relationship given by Eq. (A. 3. 9) there results 


d<f> 1 

dr 


ih / 1 - m , 2 ah 
v 2 v a? 


(A.3.15) 

Comparing boundary conditions (A.3.7) and (A.3.13) with relationships (A. 3.9) and 
(A. 3. 10), under the assumption f\ = f 2 , there result 


A = 


'p\ 


which also implies that 


|1 -M 2 \ 



(A.3.16) 

1 1 - M? 

/ 1 - Af| Q1 “ “ 2 

(A.3.17) 

1-M| 
1 -Mf 

(A.3.18) 

s related as follows 


- 

1 - Ml p2 

(A.3.19) 


(A.3.20) 



